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AliSTHAirr 

The th'hana ooherent^r.eattei' x’avlar system luiii been syntheiUsed and 
several hundred hours of echo power and line»of»sif.ht veloeity data obtained. 
The ooherent'-soatter radar utilizes a diode array previously constructed and 
components from the Urbana meteor radar. An improved receiving system per- 
mits a time resolution of one minute in the data, 

licho power from the P region shows a high degree of variability from 
day to day. lixamples of changes in power level at shorter time scales are 
also observed. Velocity data show the existence of gravity waves and 
occasionally exhibit vertical standing-wave characteristics. 
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I . INTRODUCTION 

1,1 Eavty Devetopmnte in Soatter Propagation 

The use of i*adlo waves as a means for investigating the ionosphere 
began in the early twentieth century with sounding experiments such as 
those of Appleton and Barnett [1925], A British Broadcasting Corporation 
transmitter operating at a frequency near 1 MHz led to observations of a 
i^oflocting layer at a height of roughly 85 km. Continued study of the 
ionosphere was the direct I'esult of attempts to understand and improve long- 
range communication. 

The earliest over-the-horizon communication was the familiar multiple- 
hop short-wave or HF beam transmission. Eakereley [1932] used commercial 
facsimile conununication links to investigate the propagation of wavelengths 
between Id and 50 m. Evidence in that experiment suggested a scattering 
source of "i'lnic clouds" above 100 km for these frequencies. 

World War II saw the development of higher-power, higher^frequency 
radio equipment for use in communication and radar. Reports of anomalous 
radio propagation during the war led to investigation of over-the-horizon 
propagation at microwave frequencies by the process of tropospheric ducting 
due to evaporation over the ocean. A by-product of that study was the 
discovery that the field strength beyond the radio horizon decreased more 
slowly, even in the absence of ducting, than predicted by the smooth sphere 
theory. Booker and Gordon [1950] advanced the theory that the propagation 
beyond the radio horizon is due to fluctuations in the permittivity of the 
troposphere caused by thermal instabilities. This mode of propagation is 
known as tropospheric-scatter propagation. 

Bailey et ai> [1952] predicted and subsequently discovered ionospheric- 
scatter propagation based on the theory of Booker and Gordon [1950] . A test 
path of 1245 km from Cedar Rapids, Iowa to Sterling, Virginia was used with 
transmission at 49.8 MHz, The results of Bailey et al, [1952] indicated 
that the returned signal was scattered from the D region of the ionosphere. 

Both the tropospheric- and mesospheric-scatter modes of propagation were 
widely accepted by 1955 [Norton and Weianer^ 1955] . Extensive obsei'vations 
had been made by Bailey et al. [1955] although national security interests 
limited the amount of information which could be published at the time. In 
addition, VUlara and Weiaakopf [1955] presented a theoretical analysis of 


ionosphtjric scatter, They found that turbulent mixing in the presence of a 
strong gradient electron density was responsible for the ionospheric- 
scatter propagation, The subsequent study of the niesospherlc turbulent 
regions by VHF backscattor is of interest in this paper. 

1,2 Turhulent Saatfsp Thaotnj 

Turbulent flow in a fluid occurs whenever certain conditions for sta- 
bility arc violated. In general, the motion of a fluid is governed by the 
continuity equation 

v-n = 0 (1.1) 

and the equation of motion in given by the Navicr-Stokes equation for the 
velocity field n(r, t ) ; 

{||- + Cv*A)v) + VarCVon) + 0 (1,2) 

where p, o, an<l v are the pressure, density and veloci";y in the fluid and v 
is the kinematic viscosity related to u, the coefficient of viscosity by 
V a p/p. The characteristic linear dimension with which a particular flow 
is associated is usually termed and a representative velocity is V. A 
dimensionless parameter, the Re>'nolds number i?o, associated with the above 
equations, is given by 



The magnitude of this quantity determines the nature of the fluid flow, the 
condition for turbulence being Tfc >> 1. The Reynolds number is, however, 
not the only parameter concerning tui'bulont flow. 

When a fluid is flowing in a turbulent manner eddies of different scale 
sizes are formed. The largest scale eddies are fed by the large-scale dy- 
namics of the ionosphere, namely the overall global circulation and the 
superimposed planetary waves, tides and gravity waves. The energy pumped 
into the turbulent flow must be dissipated in some manner. The large-scale 
eddies disintegrate into smallei* eddies over a short period of time and 
while doing so, dissip.ate enei'gy through viscous damping. As shown in a 
dimensional argument by VillarG and Weisskopf [1955], the amount of energy 
lost to viscous damping eventually reaches the amount of energy input into a 
given scale size of eddy. This scale is then the smallest eddy, and the en- 
ergy has been entirely dissipated. Only under certain conditions is turbu- 
lence energetically possible. Energy which is extjracted by the turbulence 
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is that which could not b6 maintninod as potential enetf)' in th**’ fluidi This 
condition for cncx’gy balance is given in terms of the Richardson nuinhor Ri\ 

m - w^vcoyaa)^ Cl. 4) 


where is the Brimt-Vaisala frequency for bouyancy oscillations, (0F^/3s3 
is the vertical shear in the mean flow, is the horiKontal velocity of the 
mean flow and a is perpendicular to the planes in which the fluid is strat- 
ifiod. Wh.cn one considers that the numerator, is a measure of the 
static stability and that (0F^/3n)^' is related to the perturbing effect of 
the wind shear tlien it is clear that a smaller Richardson number implies 
greater tendency toward instability, So, if tho Richardson number is lower 
than some critical value turbulence is onexgoticalTy possible, A critical 
value of If is often used. The above two criteria can then }>e used to eval- 
luitc tlfo nature of expected turbul.wce in the atmosphere, 

qmR Bowhi‘t7- [197S] typical values for the Rishardson and 
Reynolds numbers to examine the occurrence of turbulence, Their results 
show that one would not expect tho atmosphere to Iks turbulont as a whole 
below tho turbopausc, Instead, the Re>'nolds criterion is virtually always 
satisfied but the Richardson cri-crion can be satisfied only over regions 
small with respect to scale height. Hence, the turbulence is expected to bo 
intermittent j namely, it does not occur at all times and licights. 

Finally, the scattering induced by the turbulence must bo considered. 
The following analysis was considered by Vi.ltavs and Uoisnkopf [1955] and 
others, and was discussed most recently by Baottogi and Bowhdtt [1976a]. 
hben diffusion does not occur, the continuity equation for floctron density 
is given by; 


•dt 


^ <7 




(1.5) 


where (7 is a production term, h a loss term, and V*(i7y) a transport term. 
Production and loss mechanisms are such that g ~ L X 0, Because of the col- 
lisions of electrons with neutrals the velocity in (1.5) is taken to be that 
of the neutrals which arc assumed to form on incompressible fluid. Under 
tlicse conditions the velocity can be factored and the transport term then is 
In addition, W is generally vertical causing Z 0 since the 

vortical velocities are generally very small. It should be evident that any 
small deviation u in the velocity field will produce a ti'VN transport term 
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and hence a fUictuation ;r electron density. In this manner the turbulence 
couples to the electron density, 

An expression for the scattering cross section can bo found using 
Booker scattering theory, tot t l>e a choractoristic length scale associated 
with wave number k « charoctoristlc velocity and Gharactcriscic time 
u,, Then the fluctuations in electron density at scale t arc given by 

« (V'“ (1,6) 

where N' ^ the gradient of electron density at scales larger than t is assiuncd 
Known. Following the analysis of Raotogi and ’BaoihiVl [1976a] the length t 
is assumed small and the turbulence is assumed to bo isotropic and homoge- 
neous at this scale, Using (1.6) abovc> several results from the theory of 
statistical turbulence, and a result due to Booker the radar cross section 

cr « EOi^) rOi^) (1.7) 

is obtained where is the standard electron radius and K(k) is the energy 
spectrum of the turbulence, The above result is then used to relate the 
measured signals to the turbulence by which they arc produced . 

Turbulence is not the only cause of perturbations in the electron den- 
sity which give rise to backscattor, As discussed earlier, for small scales 
the turbulence is damped by the viscosity of the fluid, However, scattering 
from thermal fluctuations can occur producing what is now commonly called 
Thomson or incoherent scattering. A review of Thomson scatter and its use 
as an experimental tool has been written by Evane [1969] . The study of 
mesospheric turbulence induced coherent backsc.atter is an outgrowth of 
Thomson -scatter studies of the E and F regions of the ionosphere. The 
development of Thomson -scatter radars will therefore bo discussed in the 
follovdng sections. 

I,v3 Finot Indimtion of XonoopheHa Baokeoatter 

The continued development of radar tectoolcgy throughout the 1950 's 
prompted Crovdon [1958] to suggest that the study of the ionosphere by radio 
waves scattered incoherently from free electrons was possible with state of 
the art equipment. Bowlea [1958] verified the existence of this incoherent 
jicatter using a 4 -MW pulse power 40.92-^^llz transmitter in conjunction with a 
large aperture antenna at Long Branch, Illinois. Although the data collected 
at that time were very noisy the existence of a scattering region at 85 km 


'•sr 
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is quite visible in photographs token from A-scopc displays of the receiver 
output, This echo was interpreted as due to turbulence-induced scatter 
rather than incoherent scatter by virtue of its similarity to ionospherie- 
senttor propagation results in both scattering cross section per unit volume 
and fading rate [Bowlm, 1964; Blniv ct at,, 1961] , 

Additional observations by BowIob [19613 showed that the broadening of 
tlio received spcctmm for altitudes above 100 km was not ns great as pre- 
dicted by (kfi'don [lOshS] for completely incoJiorcnt scattering, Theoretical 
work such as that of PoughoPtif and Fax^lcy [19603 Indicated that the influ- 
ence of the ions on the scattering is important when the radar wavolongtli 
is larger than the Debye shielding length, The spatial variation of the 
electron density at scales greater than the Debye length is constrained to 
be that of the ions by virtue of local charge neutrality, The observed 
Doppler shift is therefore associated with the motion of the ions, The 
partial coherence in the roturned signal caused by the ion-oloctron Incei'- 
action led to the name Thoms on -scat ter, Several of the radars discussed 
below were designed without complete understanding of the Thomsoin-scattor 
spectrum effect. 

1,4 Development of Thomaon-Seattev Eadavo 

By the early 1960's several radar facilities designed for ionospheric 
study had been constructed, The radar proposed by Gordon [1958] was con- 
structod at Arecibo, Puerto Rico on the basis of fully incoherent scattering. 
The expected wide spectrum of retuCTcd signal implied low energy per unit 
bandwidth and therefore a large aperture antenna. A spherical reflector 
‘jurface of wire mesh ,^05 m in diainoter with a radius of curvature of 245 m 
is set in a limestone sink-hole, A line source is suspended above the re- 
flector and can be moved to point t(ie l>eajn 20° off the zenith in any cliroc- 
tion. A 430-Mllz 2-MW pulse transmitter is used for the Thomson-sentter 
studies. Further details of the system can be found in Gordon and LaLmde 
[1961], 

The Millstone Mill Observatory at West ford, Massachusetts originally 
employed a 440-Mll,'i tracking radar using a 2,5-MlV pulse transmitter and a 
steorable 2S-m diameter paraboloid antenna. Ionospheric -scatter measurements 
were made by Pineo ct al. [I960] using the original system, Later a verti- 
cally pointed 70-m paraboloid antenna was used with the 440-lDlz transmitter 


6 


for additional measurements while a 1295-MHz transmitter was added to drive 
the original antenna [Evane and Loemnthat, 1964]. 

The scatter radar built at Jicamarca, Peru was based in part on a sug- 
gestion by Bowleo [1961] that ion gyroresonance sidebands could bo observed 
by a. radar pointed perpendicular to the magnetic field in the ionosphere. 
Hence, the location in Peru near the geomagnetic equator was chosen for a 
49.92-MHz 4-MW pulse transmitter with a 290-m-by-290-m crossed -dipole array. 
The transmitter and antenna are designed so that the two orthogonal sets of 
dipoles can be operated independently and simultaneously to produce various 
polarizations. The modules of the crossed arrays can be fed using different 
coaxial cable lengths to produce beam steering. The facility at Jicamarca 
is described in detail by Raatogi and Bowhi'lt [1975] . 

Additional radar facilities include the French bistatic radar built in 
1962 at St. Santin and later modified to a quadristatic configuration [Bauer, 
et al.i 1974]. The transmitter operates in a continuous wave mode with 150 
kw of power at 935 MHz. The Stanford Research Institute radar operating at 
1300 MHz was originally located at Menlo Park, California [Leaddbrandt 1967] 
and later relocated to Chatanika, Alaska for auroral -zone Thomson -scatter 
studies [Leaddbyxmd et aX., 1972]. The British also operated a radar facility 
at Malvern, England on 400 MHz [Wllliami and Taylor , 1974] although this 
facility is no long'’'!' in use. 

The various radars mentioned above were used during the early 1960 's 
primarily in the study of the E and F regions of the ionosphere [Greenhou 
et al. t 1963; Evane^ 1967; Evane^ 1969]. The Thomson-scatter method was used 
to obtain electron-density profiles, electron and ion kinetic temperatures, 
and to identify the major ionic constituents by comparison with theoretical 
curves. Partial-reflection techniques were used extensively in the D region 
[e.g. Gregory, 1961]. 

1 . 5 Mesospheric Radar Studies 

The Thomson-scatter radar facilities discussed above were used for D- 
region studies beginning in the mid 1960's. LaLonde [1966] reported D- 
region 430-MHz echoes at Arecibo. VHF echoes were observed at Jicamarca by 
Flock and Balsley [1967] who found the predominant feature to be an echoing 
region near 75 km. Between 90 and 110 km the equatorial electrojet dominated 
the retLK'i'i.s. loannidis and Farley [1974] and Amisted et al, [1972] re- 
ported echoes in the region below 80 km. A coherent -scatter model was 
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required in the P region at UIIF since th6 electron density predicted by 
Thomson-scatter theory to produce the observed echo power would have been 
too largo, 

The first inonsuremonts of spectrum or autocorrelation measurements in 
the mesosphere were reported by and (hdllmi [1974’J at Jicaiiuirca, A 

velocity measurement along the line of sight of the radar beam was obtained 
from the autocorrelation function. It, was concluded that shoi^t -period fluc- 
tuations in measured wind velocities <A?ero duo to gravity-wave fluctuations 
and that neutral atmosphere turhulovice produced by wind shears was rospon- 
siblc for the observed echo power. In this investigation the signal-to- 
noise ratio of the data was improved by the process of coherent integration. 
The correlation time of the backscattored signals was longer than the inter- 
pulse period of the radar so that successive samples for a given altitude 
cotild be summed together with the signal tending to add and the noise tondi" 
to cancel, 

An o.'ctension of the above work was performed by Eaatogi and WoQihnan 
fl974] av .Ticamarca. In that investigation it was shown that the echoes 
throughout the 15 to 85 km region were duo not to Thomson scattering, i,e. 
the thermal motion of the electrons, hut rather wore turbulence induced. It 
was shown tluit the mean echo power for the received signals was 4 dh above 
that o,xpcctcd from Thomson-scattorlng theory. The authors then suggest that 
the intermiittent nature of the echoes indicates the turbulent nature of this 
region. Indeed, tlio echo power often varied by s 20 dB over a one or two 
minute period iiulicnting a nonthormal mechanism was responsible. The same 
conclusion was reached in a recent paper by Raotogi md BouMdl [1976b] . The 
large fluctmitions in the echo power arc directly related to fluctuations in 
the rate of energy dissipation per unit mass duo to turbulence. The studies 
described in this paragraph confirmed the turbulence induced origin of scat- 
tering in the mesosphoro and resulted in the construction of several rad,nrs 
for tl\o study of d>aminics in the mesosphere, stratosj)hero, and troposphere, 
(abbreviated M, S, and T) . These M$T radars utilise the coherent nature of 
the returns ami the Doppler shift produced by the motion of the turbulent 
scatcering region along the lino of sight of the radar beam. 

1 . b Mcf'>OBphotd\'j Ptmtoaplmnx't Ti'oi\mplwvia (M8T) Radm'B 

Several of tl>e radar facilities Initially used for Thomson scatter )\ave 
iKsen Giiiployod to collect coherent-scattor data from altitudes below the F 
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region. Woodimni (vid GuHtm [1974] mensurod winds and turbulence in the 
.stratosphere and mesosphere at Jicamarca. The Doppler shift of the returned 
signal indicated the motion of the scattering layer along the lino of sight 
of t!ic radar beam. Data collected at that time were limited to one height 
at a time. With the installation of a now computer, simultaneous multi height 
mesospheric obsej'vations were made and Goodman, 19771. Recently, 

waves in the lower stratosphere were observed by Rllater at ad. [19781 at 
Jicamarca. 

The 4.^0 Mils radar at Arecibo was employed by Aeo at- at. [19771 for the 
mea.surement of middle atmospheric, dynamics. As discussed by Wooiinan and 
(tuitlcn [1974] the echo intensity for coherent scatter depends on the energy 
spectrum of the turbulence. Only the spatial Fourier component of the re- 
fractivc-inde.'c fluctuations corresponding to tlie probing wavelength contri- 
butes to the scattering, Raatogt and. BowhilZ [197()b] predicted that tl»e 
wave number of the Arecibo transmitter exceeds, for mesospheric heights, 
tlic critical w.ave number above which tl»e turbulent energy spectnim de- 
creases ex|X)nontially, Aao et at, [1977] verified the prediction of Raotogi 
(md Boidhitl [19761)] but showed the existence in the stratosphere of turbu- 
lence .at the probing wave number. Additional UllF measurements at the tropo- 
pauso have l)cen made by Bats7.oy and- Fa3d,Gy [1976] using the 1290 MHz 
Chatanika, Alaska radar. 

During the recent past several new MST facilities h.ave l)een designed 
and constructed. Considerations for the design and use of MST j:adars .are 
discussed by Badalmj [l9'/8rv,b]. The SOUSY r.adar near bindnu. Federal 
Republic of Germany [Osochowahj at at., 1976] and the Sunset radar in 
Colorado [<7reen oi- at.^ 1975a, b; li^ax’noak ot al.t 1978] have obtained velocity 
measurements throughout the troposphere and stratosphere. A portable radar 
at J)okor Fl.at, Alaska was also used in tropospheric measuroment.s {Raktund ot 
di., 1977}. A permanent MST radar installation is underway at Poker Flat 
based on a prototy])e radar which is located at J’latteville, Colorado \Eektimd. 
of at., 1979]. First result.s from limited operation at Poker Flat have been 
roported [Bals'loy et at., 1980] . The SOUSY, Sunset, Platteville, and Poker 
Flat radars are all VllF radars as is the Urbana radar which is discussed 
below. 

1.7 ObJectwoB mid Saope of This Study 

The three prlncip.'il objectives of this study are: (1) to describe the 
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(lovciopinent of the Urbana coherent "scatter radar system, (2) to obtain echo 
power and radial velocity measurements with a time resolution which allows 
the observation of wave phenomena in the mesosphere, and (3) to begin assem- 
bling a data base fi:om which scientifir information can Ixs obtained. An 
outline of this study is given below. 

Cha])ter 2 dcscril) 0 s the development of the Urbana radar. The various 
hardware subsystems are discussed and the operations performed at oaci> stage 
of data processing are descrilHsd. 

Observations characteristic of the first data collected from the Urbana 
radar are shown in Chapter 3. Specific features of the data are illustrated 
and discussed in detail. 

Conclusions and suggestions for future work appear in Chapter 4. 
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2. UIVBANA COIIliRIiNT-SCA'rrriR RADAR SYSTEM 

2.1 Dtrmtopmmt of the Urbcma /fatto* 

III 1956 the National Bnroau of Standards oporotod the Long Branch Radio 
Propagation Transmitting Station Wfl near Havana, Illinois. Among other ex- 
periments conducted at this station a 40.92-MHz pulsed radar with a flat four- 
acre dipole array was utilized l>y Boidea [1958] to verify the existence of 
incoherent scatter from the mesosphere , The Havana radar was subsequently 
used by the Smithsonian Astropliysical Observatory and was made available to 
the University of Illinois in 1971. A multistatic meteor radar was con- 
structed at the Urbana field station using the d0.92-MHz transmitter [AVfotU#!, 
1974; Ueos md (hltcr, 1976]. 

Because of interest in coherent scatter fi’om the ionosphere, elements 
from the dipole array at Havana were used in a vertically pointing dipole 
array constructed at Urbana. The work on the antenna and associatetl feed 
system was completed in 1976 [A/ '/man and Bowhiltj 1976]. Collection of 
coherent -scatter data then became a problem of modific.'ition of the meteor- 
radar system. 

i\ K. Rastogi made the first observations of coherent scatter using the 
Urbana radar during July, 1976. While these verified the potential operation 
of the radar for coherent scatter, the data output was limited to oscilloscope 
A-traccs and chart recorder outputs from a boxcar integrator. In November, 
1977 the first attempt to use a computer for coherent integration and complex 
autocorrelation of real-time data was made. However, bec.ause of poor per- 
formance in various sections of the meteor -radar system being used in the 
scatter system, notably the receiver, no significant data were collected. 

After testing and repair of various components a lower noise figure 
receiving system was obtained and preliminary correlation data were collected, 
in Mvrcii, 1978. Piu'ther improvements .in the hardware brought the system to 
its present condition. Much of the hardware for tlie coherent-scatter radar 
is sliarcd with the meteor radar. Detailed descriptions of those subsystems 
arc therefore not includod in the following, but important operating para- 
liioters are repeated. 

2.2 iuuian Hanikkxi'C 

2.2.1 frmimtittcv. The transmitter utilized at Urbana for both the 
meteor radar and colierent -scatter radar was built in 1958 by Continental 
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Rlectronics Manufacturing Company of Dallas, Texas. Designed as a prototype 
for the Distant Early Warning (DEW) line of radars, the transmitter has a 
<)-MW peak pulse output power rating with nominal average output power of 
20 kW. The radar operates at 40,92 MHz with pulse widths ranging from 3 ps 
to 100 ps. Presently the system is operated with a 20-ps pulse at n pulse 
repetition froquency of 400 Hz . 

The final amplifier stage consists of four water-cooled triodo tubes 
(Machlott ML-5682) in grounded-grid configuration. The tubes are mounted in 
pressurized cylindrical cavities to minimize arcing. The cavities are reso- 
nated at the input and output of the tubes. Water cooling is provided to 
maintain operating temperature near 100® F. A cavity pair is shown in Figure 

O 1 
• .1 • 

The transmitter final stage is actually two pairs of tubes with gach 
pair driving one side of a balanced coaxial line. The driver amplifier, 
delivering 0,6 peak pulse power, is a single ML-5682 triode configured in 
the same manner ns the output tubes. Power from the driver is split first 
to obtain the 180® phase shift between pairs of final tubes, then again to 
properly phase the drive to each tube in the pair. 

The intermediate power amplifier is a 4CX5000 tube which operates at 
reduced plate voltage with respect to the ML-5682 tuixjs. Plate dissipation 
of S KW necessitates a pressurized chamber with water-cooled walls. The 
4CX5000 is in turn driven by a Continental Electronics 814B VHP transmitter. 
Using a pair of 4CX1000 A's connected in parallel the 814B provides drive 
pulses with a peak power of 3 kW. Various lower power stages drive the 814 B 
from the gated RF pulse provided by the radar director. 

The pulse modulator produces amplified, shaped pulses which are trans- 
former coupled into a final switching stage. Three parallel ML-5682 tubes 
switch current in the pulse modulation transformer, producing the required 
high voltage pulse on the secondary. Partial control of the transmitter 
output power is obtained by controlling the step-up ratio of the transformer 
The final switching stage is shown in Figure 2.2. A more detailed descrip- 
tion of the transmitter can be found in Heae and GeZleVt [1976] . 

2.2.2 Antenna, and tvmiamit/veaeive euH-tah. In conjunction with a 
powerful transmitter coherent scatter radar requires the use of a high-gain 
narrow-beam antenna. The characteristics of the Urbana dipole array are 
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•l().l)2-N0lz transmitter. In the foreground of the figure 
are the control panel, the power meters and the circuit- 
breaker panel. The lar^e metallic cylinder in the center 
of the picture is a pressurized, tuned cavity housing, the 
driver tube. One of the four output cavities is seen in 
the backj>round on the left. The equipment located im- 
mediately in front of the output cavity in the fiRure is 
the transmitter for the lirhana laser radar. 
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r<nal switching stage of pulse modulator. The 
three large tubes seen here are Ml.-SfoS.. t nodes. 
Ml.-5b82 tubes are also used in the driver an 
power amplifier cavities. 
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given in Tnblc 2.1 [A'tV*ia>\ <xnd Bowhil'l, 1976], Tho array is organised into 
groups of twenty-eight pairs of half-wave dipole elements called cell.s, Six 
cells arc connected in parallel to form a group. Tho 1008-element array con- 
sists of six groups. The feed system for the antenna employs open-air trans- 
mission line transformers for iiuatching. The feed line from transmitter to 
antenna is 450 ft, of balanced coaxial line built from aluminum construction 
tube* The inner conductor of three inch diameter and oute'i^' conductor of seven- 
inch diameter arc separated at seven- foot intervals by teflon spacers. The 
antenna is shown in Figure 2.3 while the feed line is pictured in Figure 2.4. 

The antenna is connected to the transmitter via a transmit/recoive (T/R) 
switch to allow monostatic operation. Four gas-filled tul)es are employed 
forming a txansmit/roceive switch for each side of the balanced coaxial lino, 
The transmitter voltage causes breakdown in tho tubes and the resulting low 
impedance is transformed so as to isolate the receiver during tho transmit 
pulse and the transmitter while receiving. Recent tests measured the recovery 
time of the switch at 400 yis when operating witli l-MW peak pulse iwwer. De- 
tails of the T/R switch, antenna, and feed system design and construction are 
given, by Attwan and Bowhitl [1976]. The transmit/receive switch slicd is shown 
in Figure 2.5, 

2,2.3 Eadav tUreatov. The radar director at Urhana consists of two 
separate units; one for radio frequency (RF) s>aithesis and one for timing and 
pulse generation. The RF section contains crystal oscillators at the trans- 
mitter .frequency of 40,92 Mllz, and the receiver local oscillator frequency of 
35.42 Mllz. Mixing these two frequencies and phase shifting by :'’45'' yields 
quadrature reference signals at the intermediate frequency of 5.5 Milt, 

The timing section of the radar director is drivcMi by an external master 
clock of 100 kl.lt derived from a. 5-MMt reference, livery pulse intcn^^l is a 
multiple of tlie base lO-gs period. The transmitter, analog-to-digital con- 
verter (ADC), blanker, and the oscilloscope displnys are all driven by 
various pulses from the radar director, A detailed de.scription of the asso- 
ciated circuitry is given by Ile&a and Gettop [1976] . Timing information for 
tlie radar dire t03r is given in Appendix H, 

'2,2,4 Receiving syotenu The receiving system for the Ui’bana. coherent - 
scatter radar shares a. receiver and quadrature phase detector with the meteor- 
radar system. To protect the receiving system and lower the noise figure, a 
blanker/preamplifier unit which can he locatetl near the T/R switch has been 


Table 2.1 

Antenna parameters for the Urbane coherent-scatter radar . 


Aperture illumination efficiency 
Antenna efficiency (p^) 

Radiation efficiency (p^ = hw/h.,*) 
Physical aperture (A ) 

Effective aperture (p A ) 

(X ^2 

Directivity « 4ir«4^p^/A 0 
Power gain {g^ » p^^) 


0.69, -1.6 dB 
0.17, -7.6 dB 
0.25, -6.0 dB 
11000 
1870 m^ 

1800, 33 dB 
450, 27 dB 








Fcedlinc for the coherent-scatter system. Twin 
seven- inch coaxial lines run approximately 450 feet 
from the switch shod to the antenna. The meteor- 
radur antenna can be seen at the right of the figure. 
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Figure 2.5 bottom view of switching shed. The output to the 
coherent-scatter antenna is at the right center of 
the figure while the output to the met .'-or- radar 
antenna is at the upper left. The metallic bands 
in the center of the foreground delineate one of 
the four tubes used in the transmit/receive switch. 


constructed, The blanker utilizes PIN diodes to provide low Insertion loss 
and higiv power handling eapability, A low noise preamplifier immediately 
follows the blanker and essentially determines the noise figure of the sys« 
tern. The blanker/preomplifier unit is located roughly 60 feet from the 
receiver in order to minimize the degradation in noise llgurc due to coaxial 
cable between the T/R switch and receiver. Information on the blanker is 
contained in Appendix I, 

The receiver is single conversion j«lth a bandwidth of 230 kHz centered 
around d0,92 MHz, A toggle-switch attenuator is located between the mixer 
and the IF section to provide control of .signal levels at the $,$ kHz IF 
frequency, The receiver output is amplified and applied to two four-quadrant 
multiplier chips which also receive the appropriate quadrature phase refer- 
ences at 5.5 MHz generated by the radar director. The quadrature-dctoctod 
signals are filtorud to reduce the bandwidtli to 75 KHz, An additional out- 
board filter is added for coherent scatter to bring the bandwidtli down to 
40 kHz, A detailed description of the receiver and phase detector is given 
by //caa and Oeltot\ 1976] , 

2,3 Data D^oaenaing Hardware 

Analog-to-dlgltal conversion for the Urbona radar Is performed by a 
Hewlett Packard 5610 converter, The .5610 is a 10-bit A/D converter operating 
at 100 KHz with an accuracy of A'ihSB, A 16-channel multiplexor is also pro- 
vided. Control of the converter is performed by an outboard sequencing 
circuit and several mode switches on a custom built interface. 

The majority of the data processing for the Urbana radar system is accom 
plished by the Digital Equipment Corporation PDP-15 computer located at the 
field station. The PDP-15 contains 52K of core memory, an extended arithmeti 
clement for hardware multiply and divide, a real-time clock and a higli-spood 
data channel for input/output . In addition, four DriCtapc units provide 
medium speed bulk storage while four fixed-head disks are u.sed for high- 
speed storage and system software, Data tTansfer to other computers can 
lx; accomplished using the high-speed paper tape rcader/puncli. A General 
Electric Tcrminet 1200 teletype producing hard copy at 120 characters per 
second and an Infoton cathode-ray terminal operating at 9600 baud provide 
the operator interface. 

The plotting for the coherent-scatter experiment is performed on u 
Hewlett Packard 9830A desktop computer. The 9830A is a small, slow general- 
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purpose machine which utilizes a BASIC interpreter in ROM. Additiona'!. ROM 
chips for string variables, matrix operations, paper tape reader control and 
plotter control extend the capabilities of the computer. The 9830A at the 
Aeronomy Laboratory contains 16K of RAM and a single digital tape cassette 
unit which is an integral part of the computer. Digital cassettes are used 
for both program and data- file storage. 

2 . 4 Softmre 

2,4.1 Theory. The purpose of the system software is to produce echo 
power and velocity measurements at the desired sample heights. The power 
and velocity information can be obtained from either time domain or frequency 
domain information by finding respectively the autocorrelation function or 
the power spectrum. The frequency spectrum approach utilizes narrowband 
filters to obtain the energy at various frequencies. The output of the 
filters is then sampled and digitized. In the time domain approach used 
at Urbana the phase detector outputs are sampled and autocorrelated digitally^ 

The phase detector output consists of both noise and signal components, 
I’Urthermoro the noise and signal are uncorrelated so that the autocorrelation 
of the sum of signal and noise is the sum of the autocorrelations of each 
component taken separately. The signal energy is concentrated about the 
Doppler frequency of the returned signal while the noiso has a band- limited 
spectrum. Thus the signal is correlated for times much longer than the 
correlation time of the noise and as Goodman and Guillen [1974] showed the 
signal is correlated for times longer than the interpulse period while the 
noise is not. The process of coherent integration where the samples taken on 
successive pulses are added can thereby improve the signal -to-noise ratio of 
the resulting sum. The coherent integration period is chosen as a compromise 
between improvement in signal -to-noise ratio and time resolution in the auto- 
correlation process which follows. 

The power in the returned signal can be calculated from the zeroth lag 
of the signal autocorrelation. If B(r) is the autocorrelation of a process 
37(1) then considering the definition 

B(r) = S(z'Ci + t;) (2.1) 

one has with t = 0 

R(a) - E{x(t) = E(lx(t)j^} (2.2) 

which is an expression for the power. However, as pointed out above the 
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noise will contrihuto primarily at the zeroth lag of the total autocorrela* 
tloti. If the noise power is assumed constant over the sampled altitudes and 
slowly varying then changes in power calculated from the zeroth lag of the 
total autocorrelation function will l)c due to changes in signal power. 
Relative fluctuations in power over altitude and time ni-e therefore moasurod. 

The motion of the scattering regions can he measured using the Doppler 
shift of the returned signal, If is the r.adial velocity of the scattering 
region and ,f| is the Doppler frequency then 
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whore X !=> 7.53 meters is the radar wavelength. Hence 
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C2.33 


C2.43 


where is the radian Doppler frequency, and <}t is the phase of the autocor' 


relation function. To estimate ^ one recalls th.at (fCo) « 0, i.e. the auto- 
correlation is real at the zeroth lag by definition. So one has flie following: 

^ I 

*1 e 0 


1 


(2.5) 


A wciglvtcd average of the values calculated at lags one to three is used to 
olvtain the velocity measurement. 

2.f|,2 Data coVle&tion. The first state of d.ata processing for the 
IJrbnna rnd.ar consists of data colleetion and reduction to autocorj'clation 
functions in real time. Tl\ere are four distinct processes which occur in 
real time: input, colicrent integration, correlation and averaging. These 

four tasks are interlaced by the software in such a way tl>at the latter three 
arc interrupted whenever the analog to digital convertor is sampling. The 
characteristics of the data collection process are discussed below. 

The aTtitude region of interest for the coherent-scatter experiment is 
from 60 to 90 km, or a 50 km inteiafal . The altitude interval detcjininos the 
number of samples In the following manner: the distance to a scattcrer is 
detennined by iiioasuring tlie time taken by the radar pulse to travel to the 
scatterer and rcturji. One has 
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where R As the range, a the speed of light and t the time for the wave to 
travel up and back. A time difference of At « *2 “ Iwtween sample pulses 
results In a range resolution Ai? given by 

^^(*2 - j ) 


m ^ 


(2.7) 


The 10-pscc conversion time of the analog-to-digital converter therefore 
corresponds to a range resolution of 1.5 km and 20 samples ai’e required to 
cover the .50 km region of interest. 

The cosine and sine channels of the phase detector are sampled on alter- 
nate radar pulses requiring a pair of pulses to obtain a complex sample at 
each of the 20 heights. With the pulse repetition frequency of the radar at 
400 lla a pulse pair requires 5 msec. Twnty-five such pairs of pulses fill 
an input buffer after 1/8 second. The input is double-buffered so that data 
collection can continue on an intcjrrupt basis while processing occurs. 

The coherent-integration process adds the corresponding complex values 
at cacli lieight reducing the twenty- five complex sample sets to a single set 
of twenty complex values. Note that the twenty- five additions can add no 
more tlian five bits to the ten-bit words from the a ralog -to -digital converter 
so that fifteen-bit wojxls are produced. Single-precision addition which 
saves processor time is therefore possible on the 18-bit PDP-lS. The process 
of colierent integration As illustrated in Figure 2.6. 

The data from the most recent 1/8-second coherent integration interval 
is correlated with data from the previous intervals. Consider a stack of 
data sets, each set containing twenty coherently integrated complex samples. 
After each l/8-second period the latest set of samples is pushed onto the 
top of the stack into the pi'esent time slot while the oldest data arc lost 
at the bottom of the stack. Thus the complete complex autocorrelation for 
all lags and for each of the twenty heights iiaist be completed during the 
1/8-second colierent integration time. The speed of the PDP-15 allows a 
multibit-multibit complex autocorrelation to be performed out to lag 12, ot- 
to 1.5 seconds. The output of the complex autocorrelation function is a 
real value at lag zero, and real and imaginary values for lags 1 to 12. 

Twenty- five numbers arc thus produced for each of the twenty heights every 1/8 
sccoiul. The data are now double-length words as a result of multiplication 
durin g autocorrelation . 

Tlic autocoiTClation functions are then converted to floating point and 
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aro averaged for one minute by adding. Four of the one-m,inute arrays formed 
by averaging are stacked to produce an output buffer which becomes a single 
record on disk. Three records form a file, and ten files fill the disk. 

There are three disks available so that after six hours the data collection 
terminates and the averaged autocorrelation functions are either processed or 
stored on DECtape for later processing. Data-collection programs are listed 
in Appendix III. 

2.4.3 Vast ppGoeaeing and data tvanefer. Upon completion of data col- 
lection the correlation functions stored on DECtape or disk are processed 

to obtain useful scientific information. At the present time values for the 
returned echo power and the velocity of the scattering region are calculated 
at each sample height for every minute . The output of these programs is to 
paper tape so that data may be transferred to the HP-9830 for plotting. 

The power and velocity are calculated from the autocorrelation function 
as described previously. The logarithm of the power is punched onto paper 
tape to increase the dynamic range of the data which can be transferred. The 
velocity program tests the autocorrelation fvmction for noise characteristics 
by comparing the amplitude of the autocorrelation at lags 1 through 3 to the 
real part of the autocorrelation at lag 0. If at a given lag the ratio is 
too small then that lag is not used in the velocity calculation. If all 
three lags fail the test a value of zero is punched onto the velocity t.ape 
as a signal to the plotter software that no velocity was calcul.ated at that 
point. The post-processing programs are listed in Appendix IV. 

2.4.4 Blotter software. The plotting of data from the coherent- scatter 

radar involves two stages: (1) transfer of data from paper tape to cassette 

tape, and (2) plotting of data from the cassette tapes. To store the data 

on cassette a paper tape containing two hours of power or velocity data are 
read into the HP-9830. Additional information which identifies the data are 
tyijed into the 9830. A cassette file is then formed from the data and iden- 
tifying information. The use of cassette storage allows multiple use of the 
data file without the necessity of rereading the paper tape. 

Three types of plots are made for each data set of two-hour duration. 

The velocity plot is one of velocity vs time at a givf'n altitude. When the 
data exceeds an operator-controlled limit value or when the value zero occurs 
the data point is not plotted. Recall that zero was the signal value gener- 
ated by the PDP-15 to indicate that no velocity value was calculated because 
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of the noise properties of the autocorrelation function. Two types of power 
plots are produced, the first being power vs altitude at a given time or 
power profiles. The power profiles show the variation in altitude of the 
scattered power but the temporal variation is best observed on the second 
plot, the power vs time at a given altitude. Both power plots employ a 
hiding routine to improve the readability by preventing crossings of the 
traces. Sharp increases in the power level which last for one minute are 
almost exclusively due to meteor echoes. These bursts of power tend to hide 
the scattered power Values on the power profiles. A limit routine \4tich 
clips one minute spikes to an operator controlled value helps minimize the 
meteor contamination of the power profile. Examples of these plots are shown 
l>elow in Chapter 3. The plotting routines are given in Appendix V. 
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3. MESOSPHERIC COIE RENT- SCATTER OBSERVATIONS AT URBANA 

3.1 Obeewation 'Pvogvam and Introduatovy Rmanke 

The Urbana coherent-scatter radar was first operational on a day-to-day 
basis in April, 1978. Extensive data collection occurred in that month often 
on a sunrise to sunset basis. The collection of data continues at present so 
that a data base can be accumulated. As of January, 1979 242 hours of coher- 
ent-scatter data from the mesosphere has been collected on 34 days. 

The figures shown in the sections which follow are not above average in 
their measurement quality. No smoothing is employed to obtain idiat appear to 
be continuous curves: rather the time resolution of the system, one minute 

as mentioned earlier, is sufficient to yield the plot quality obtained. The 
figures, however, are chosen to illustrate in a consise manner features of 
the observed data . Data from Urbana of a similar nature has been shown in 
Milter et at. [1978] . 

3.2 Echo Power Data 

Power data from April 21, 1978 is shown in Figure 3.1. The presentation 
in this figure has been referred to earlier as power vs altitude at a fixed 
time, or power profiles. The areas of no returned signal power, generally 
below 66 km but not limited to this region, appear as essentially flat, 
evenly spaced lines. The bumps in the profile therefore exhibit power levels 
above the noise. In Figure 3.1 the power from the scattering regions is 
roughly 3 to 5 dB above the noise with the possible exception of the scat- 
tering region at 87 km at the end of the two-hour period. Activity such as 
is shown would be considered characteristic of a quiet day. 

Note further the nature of the scattering regions during the two-hour 
period. A region centered at 69 km is steady for the entire period. A 
second layer near 73.5 km seems to drift slowly downward perhaps merging 
with the layer below by 1030 CST. Finally, a third distinct echoing layer 
appears at 79.5 km near 940 CST. The scattering regions in Figure 3.1 are 
well defined, change slowly or not at all in altitude, and can appear and 
disappear abruptly. 

A more active day, April 11, 1978 is shown in Figure 3.2. The echo 
power is often 10 dB above the noise at altitudes above 69 km and scattered 
power is observed even at the lowest altitudes. The scattering regions 
appear more mobile particularly at the lower altitudes. Furthermore, it is 
not clear when the lowest layer near 63 km has disappeared. At 915 CST the 
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Vertical profile of scattered power at Urbana beginning at 854 CST 
on April 21, 1978. 





P.RX= 12.03 BflSE= 7.23 



Vertical profile of scattered power at Urbana beginning at 800 CST 









layer begins to merge with the layer above at 72 km but at 930 CST briefly 
reappears as a very weak scattorer only to be gone again by 945 CvST, The 
high scattered power values obtained in Figure 3,2 are due to a solar flare 
which occurred shortly before the data set shown, at 740 CST. 

An unusual example of variability is shown in Figure 3.3 from May 18, 
1978. Note the slightly higher altitude of the plot exposing the region 
above 87 km alluded to in the previous examples. The highest region exhibits 
wide changes in scattered power. Furthermore, a lower power region directly 
below, at 82,5 km, seems to increase in power out of phase with the region 
above, 'file lower scattoring layers, near 72 km, do not show the degree of 
change of the higher altitude layers. 

3,3 Vatooity Data 

The wave nature of ionospheric dynamics is observed in Figures 3.4 and 
3,5. Tlie velocities shown hero correspond to the power plots of Figures 3.1 
and 3,2 respeGtively . Velocity data for April 21, 1978, Figure 3,4, indicates 
low amplitude oscillatory behavior particularly at 930 CST. Comparing the 
velocities with the power data in Figure 3,1 one observes that when the re- 
turned power is high the velocity curves are continuous. More precisely, if 
the autocorrelation function is very noisy due to low returned signal power 
then the velocity obtained from the autocorrelation function is not plotted. 
Generally the one- minute time resolution is adequate to produce smooth veloc- 
ity curves. When higher frequency waves are present with high relative am- 
plitude however, the rough velocity curve observed at 87 km is obtained. 

Note also that the mean value of the velocity for the two-hour i^criod is 
negative, particularly at lower altitudes. The non-zero value of the mean 
indicates a horizontal component of velocity corresponding to the off-vertical 
pointing direction of the radar. 

Higher amplitude velocities are observed in Figui*e 3.5. Again waves 
are readily observed, with an eight-minute period wave the predominant com- 
ponent. It should be noted that for the two examples given, the higher 
amplitude velocities correspond to the day with greater overall scattered 
I>ower. Generally the more active days tend to produce higher velocities; 
however a few counter examples have been observed. Hour- to-hour variations 
in the scattered power do not appear to affect the nature of velocities 
obtained but rather the continuity of the plot as explained above. 



Vertical profile of scattered power at Urbana beginning at 1114 GST 
on May 18, 1978. 





velocity at Urbana beginning at 854 CST on April 21, 1978 
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M ndditional fonturc of tive velocity data is illiisferated by a plot 
from April 13, 1978, Figure 3*6. Comparatively high velocities are observed 
throughout the two-hour period with oxamplos of sevoral froiiiuericies of waves 
appar*ent, Of special interest in this figure is the evidence for vertical 
standing waves. Consider 1107 CST and examine the velocity as n function of 
altitude. The velocity changes sign at this time between 81 and 79,5 km am! 
again at 69 km. The same behavior is reodily observed at 1150 and 1210 CST, 
Similar ovidence for vertical propagation of gravity waves occurs in other 
data sots. 


3.4 ObBoinmtims of an Entire Pay 

The data sot.s discussed above are two-hour segments which illustrate 
important aspects of the observations* Any physical procosses which occur 
at longer periods however require the study of observations taken throughout 
the day. For this purpose the complete set of observations for May 24, 1978 
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to Figure 3.21. ITio plots are in throe groups: 


1, the power profiles, 2. the plots of power vs time at fixed altitude, and 
5. the velocity plots, 'the data is continuous from 812 CST to 1820 GST save 
for a seven ’•adniito transmitter failure at 843 CST and oiglit minutes at 1412 
CST required to dump the data onto tape. 

Comparison of Figures 3.7 and 3.12, the two power plots beginning at 
812 CST, shows the iiwinnor in which the power plots complement each other. 

Tlie power profiles in Figure 3.7 at 900 CST and 87 km are hidden so that 
one must use Figure 3.12 to dotermino how quickly the power level decreased. 
In addition, the software iiiniting routine apparent at 925 CST in Figure 3.7 
is sot to a low value so that the nmximum roadability of the power profiles 
is obtained, Ihe same limiting routine is avail.able on the power vs time 
at fixed altitude plots but is generally not employed so that the full value 
of one-niinute peaks appear as in Figure 3.12, From the second plot one can 


judge if data have been clipped and if so the limit value is changed and a new 
plot obtained. 

Tlie data from May 24, 1978 show an active day. The velocities arc high 
in amplitude and show many different frequencies, some with periods as long 
as twenty minutes. From 1315 to 1415 GST evidence of vertical propagation 
can be observed as wave crests occur at progressively later times for diffox- 
ent altitudes. Finally, consider the scattering region between 84 and 94 kin. 



elocity at Urbana begiTining at 1016 CSX on April 13, 1978. 
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Vertical profile of scattered power at Urbana beginning at 1212 CST 
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Vertical profile of scattered power at Urbana beginning at 1421 
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Horizontal profile of scattered power at Urbana beginning at 1012 CST 
on May 24, 1978. 
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Horizontal profile of scattered power at Urbana beginning at 1212 CST 
on May 24, 1978. 
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Figure 3.18 Line-of-sight velocity at Urbana beginning at 1012 CST on May 24, 1978 













velocity at Urbana beginning at 1621 CST on May 24, 1978 





Tn Figure 3,7 the region is not so well doFined as later in the day yet still 
appears to bo shifting in height. Throughout the day the layer moves slowly 
downward then in a shorter span of time moves back to the higher altitude and 
begins descending again. The slow descending process is most obvious from 
1421 to 1621 CST. The layer appears to split before fading away at 1630 CvST, 
Again at 1750 CST the layer is descending during it's brief reappearance. 

'Ihc long term motion of the scattering regions as shown in Pigui'cs 3.7 to 
3,11 is seldom seen at the Urbana radar. 
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4. SUNWARY AND SUGGESTIONS FOR FUTURE RESEARCH 

4 . 1 Sumnary 

The principal conclusions of this study are summarized below. 

(1) Tlie Urbana col\erent -scatter radar has been synthesized from the 
meteor radar by modification of the radar director and other meteor-radar 
components. A flat dipole array is used in conjunction with a 4-MlV trans- 
mitter operating at 40.92 MHz and a receiving system connected to a minicom- 
puter to obtain echo power and velocity measurements. 

(2) Tlie sensitivity of the Urbana radar allows data collection through- 
out the mesosphere on virtually all days. A time resolution of one minute is 
obtained which produces continuous plots of velocity without smoothing. 

(3) Tlie 242 hours of data collected from the mesosphere show the varia- 
bility of scattering. Both active and quiet days are observed with echo 
power exhibiting a dynamic- range of 20 dB. Gravity waves are observed with 
evidence for vertical standing waves. 

4.2 Suggest-ions for Future Reaearah 

4,2.1 Improved range reaoZution. Tlie turbulent layers which are re- 
sponsible for coherent scatter have a vertical separation on the order of 
several kilometers , but a thickness on the order of tens to hundreds of 
meters as shown by Rastogi and Botdvill [1976b], It is therefore desirable to 
obtain a range resolution of one kilometer or less so that the scattering 
regions can be studied separately. Detailed study of the vertical structure 
of the velocity field will also require improved height resolution. 

Tlie range resolution of a pulsed radar is limited by the pulse width and 
the associated receiver bandwidth. Decreasing the pulse width improves the 
range resolution: but if the pulse repetition frequency is not also adjusted 

then the average transmitted power will decrease thereby reducing system sen- 
sitivity. Coding the radar pulse allows the transmitter to be run at maximum 
average power with a long pulse while obtaining the range resolution corre- 
sponding to a fraction of the pulse width. The phase of the transmitter is 
varied during the pulse according to a code word. The minimum time between 
phase shifts determines the range resolution. The returned signal is decoded 
by correlation of the returned signal with the transmitted code. 

Implementation of coding at Urbana will require several modifications to 
the system hardware. First, a radar director capable of producing the coded 
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low level RF for the traJismitter must be constructed. The radar director 
would also control the decoding process. Second, the increased range resolu- 
tion will require nore altitude sampling bins to cover the entire mesosphere. 
Ihe data rate will exceed the I/O capability of the PDP-15 and the present 
A/D. Furthermore, the CPU time required for decoding will not be available 
if more sample heights are used. A decoder-preprocessor external to the 
PDP-15 and under direct control of the radar director could decode the in- 
coming signal and perform coherent integration. The PDP-15 would then have 
CPU time available for computing the autocorrelation functions for the 
greater number of altitudes. 

4.2.2 Measurement of hox^usontal veloaity. A monostatic radar can 
measure the velocity of the scattering medium along the line of sight of the 
antenna beam. If the beam is not strictly vertical then the horizontal com- 
ponent of velocity in the direction of beam tilt can be measured. Presently, 
the Urbana radar beam points about 1.5 degrees from the vertical in a gener- 
ally southeasterly direction. Non-zero average values for velocity are there- 
fore interpreted as horizontal velocity toward the southeast. By shifting 
the pointing direction to due south and again to due east two orthogonal 
horizontal velocity vectors could be obtained although these velocity values 
would be from different scattering volumes. In general, any two pointing 
directions not along the same line from the vertical can be used to obtain 

the orthogonal velocity components. 

Tlie Urbana array consists of three modules along a northeast/southwest 
line. As discussed earlier the modules are presently fed in phase which pro- 
duces a pointing direction of roughly 1.5 degrees from the vertical in a 
southeasterly direction. Tlie non-vertical direction of the beam is a direct 
result of the ground slope beneath the antenna. By feeding the modules with 
different phases the beam would be steered along the northeast/ southwest line 
and two orthogonal pointing directions could be obtained. 

4.2.3 Additional data processing. Data collected at the Urbana radar 
is routinely processed to obtain plots of echo power and line-of-sight veloc- 
ity. The minute-by-minute variation of these parameters is readily observed 
but comparison at longer time intervals is nwre difficult. Obtaining one 
hour statistics from the data would facilitate long teim comparisons and aid 
in the study of possible relationships between scattered power and observed 
velocity. 
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TTie utility of a one-hour average of velocity data has been discussed 
above. In addition to obtaining the average one would like to have parameters 
which indicate the relative amplitude of the velocity waves and the dominant 
wave period. The relative amplitude of the waves from hour-to-hour can be 
observed by calculation of the second moment about the mean, the variance of 
the line-of-sight velocity. It is important to subtract the mean here because 
the mean is assumed to be due to horizontal motion while the short-term vel- 
ocity variation is the vertical component. The standard deviation, the square 
root of the variance, therefore provides a root -mean -square estimate of the 
wave amplitude. A parameter related to the dominant wave period may be cal- 
culated in either the frequency or time domain by calculating the Fourier 
transform or the autocorrelation function respectively. 

One hour statistics for the scattered power serve primarily to relate 
the power data to the velocity information. A single power profile summa- 
rizing an hour of data would therefore be adequate but should not be obtained 
by averaging. Observations of scattered power at Urbana, as illustrated 
earlier, show short bursts of high returned power which are primarily due to 
meteor echoes. An average of the power data would thus be affected by the 
returns due to meteors. Calculation of a median value of power, however, 
would reduce the effect of short duration, high power returns and produce the 
desired characteristic power profile. 
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APPENDIX I RECEIVING SYSTEM BLANKER 

The Urbana radar is a tnonostatic system and therefore requires a means 
of isolating the transmitter and receiver. The protection system consists 
of the high-power T/R switch anJ'the low-power blanker illustrated below. 
During the transmit pulse the T/R switch limits the amplitude of RF on the 
receiver system input to roughly 60 volts across 50 ohms. The protection 
afforded by the blanker is therefore necessary to prevent damage to the 
receiving system which follows, The blanker/preamplifier unit is shown in 
Figure Al.l. Schematic diagrams for the drive circuitry and for the blanker 
itself are given in Figures A1.2 and A1.3 respectively. 





Drive circuitry' for the receiving system blanker. 
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APPHNDIX II, RADAR DIRECTOR TIMING INFORMATION 
'Fhc circuitry oinployed in the radar director is shown in HeBo and Gcllor 

[1976] . However the operation of. the radar director was not described in 

that work. 'Phe definitions and timing diagram below explain the function of 

the various controls and the resultant pulse ti*ains. 

BASIC TIME PFsRIOD: The period of the main time base driving the radar direc- 

tor, A 1-MHz clock is divided down to 100 kHz to produce a 10 psoc basic 
time period. 

PRF WORD: A thumbwheel switch which determines the number of basic time 

periods between the leading edges of consecutive PRF pulses, T}ie PRF 
word switch is set to 250 for the scatter radar producing a 400 Hz 
pulse repetition frequency. 

RANGE WORD: A thumbwheel switch which determines the number of basic time 

periods between the leading edge of the PRF pulse and the following 
Range pulse. The range word should always be less than the PRF word. 

Tbe range word is set 249 for the scatter radar which results in a 
Range pulse beginning 10 psec before the PRF pulse. 

START/STOP: A thimibwheel switch which determines the number of sample pulses 

and when with respect to the PRF pulse they occur. The 4 switches on the 
left determine the start time, those on the right the stop time. The 
difference between the stop and start settings is the number of sample 
pulses. Tlie first sample pulse occurs at Start +1 basic time periods 
after the leading edge of the PRF pulse. Start should be less than Stop. 
Settings generally are 40 and 60 for the coherent- scatter radar. 

PULSE WIDTH: A set of toggle switches which determines the length of both 

the PRF and Range pulses. The pulse length is the number of basic time 
periods corresponding to toggle switches in the up position. The LSB 
is located at the left. A setting of seven or greater is used for the 
20 iiscc transmitter pulse length. 

PRF PULSE: A positive going pulse which drives the transmitter circuitry. 

RANGE PULSE: A positive going pulse which drives the RF gater, blanker, 

and other equipment. 

ECHO SAMPLE WINDOW (ESW) : A negative going pulse train which drives the 

analog to digital converter. The width of the pulses is that of the 
driving time base pulse and the separation between pulses is the basic 
time period. 


jLJUuixyjuuuuie^ 


main time base 


_r 


RANGE PULSE 


-2^ 




PRF PULSE 




■ih 


ESW 


I RANGE WORD x BASIC TIME PERIOD 


IPRF WORD X BASIC TIME PERIOD 




1 1 

j I 

I 


I (START l)x BASIC TIME PERIOD 


STOP X BASIC TIME PERIOD 


Figure A2.1 Timing diagram for the radar director. 
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APWjNDIX I II DATA COl.MiCiriON PROGRAMS 
The two programs boiow> PSCAT mul DM ore the real- time colloction 
progroms Tot* the Urbuna coherent»scattor radar, The program DM is 
required by program siiso constraints in the PDP-15, FSCAT is used to 
perform all the data col lection and to produce correlation functions of 
12 lags of 1/8 second for each of 20 sample altitudes at one-minute 
intervals, Twelve minutes of data form a file and ten files fill a disk. 
Data collection must bo intemiptod after six hours to empty the disks 
onto DRGtapo. Before running the program the file pointer must be edited 
to indicate the correct date, The edit command "1, .SIX" is used to find 
the appropriate program lines. The date are then changed in the file name 
and the file extension shown in the line that follows may also have to 
bo initialised with a new value, The modified program is compiled with 
the mCRO compiler and then ti-'ansferrod back onto DRCtape. The .DAT 
slot assignment for loading and executing the programs is as follows: 
disks 1, 2, and 3 are assigned to slots S, 6, and 7 respectively, and 
the Dl-Ctape unit containing the FSCAT and DM binary files is assigned to 
slot -4. Disks 1, 2, and 3 should be nulled before loading and execution 
begins. 


DM program 


CM 


.GLOBL CM 
.BLOCK 764z 
. ONO 


FSCAT program 


.GLOBL MAIN, CM,. da, .PA 


.IQDEV 

5,6,7 


SHAL«6O0CJ00 



CLnN»/ 00044 



CL0F«7 00004 



MAIN DZM 

SEQM 


LAC 

CM 

/SET ADDRESS IN WRITE STMT 

DAC 

BOUT+? 


DAC 

WATEt?- 


Tad 

(7640 

/POINT TO SEQ# ADDRESS 

DAC 

SEQaD^ 


I AC 



DAC 

WDAD# 

/AND SWITCH address 

.ENTbR 

5, FI Lb 


CIDCK .'IIMBR 

Z00, ADERR,7 

.timer 

0, SYNC, 6 


. IDLE 



PADS 0 


/ENTRY 

.close 

5 


LAC 

SEQ 


JMS 

.Pa 


ISZ 

EXT 


GOO . enter 

5, FILE 


CLON 



.TIMER 

0.SYNC, 

6 

.RLXIT 

paus 


SYNC 0 


/RESTART ENTRY 

LAC 

SEQ 


SNALCLL 


I AC 


/FAKE )T FOR FIRST RECORD 

LMQ 



CLA 



DIV 



3 


/DIVIDE BY 3 

SZA 


/TEST REMAINDER 

JMP 

LEFT 


L .CLOSE 

5 


ISZ 

EXT 


LAC 

SEQ 

/TEST TO CHANGE DISKS 
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FSCAT program (cont.) 


sad (36 

JMP CLOP 

SAD (74 

JMP CLOP 

Sad ( I 3z 

JMP CLOG 

3 .ENTfR 5, FILE 

LEFT JMS A DC SET 

JMP .+3 
.DSA ONE 
.DSA RETURN 

. idle 
samp 0 

ONE I 

TIMERR 0 

N 1 750 

STOP 0 

RETURN RETURN+700001 
0 

LAC BUFFI 
and (740000 
SZA 

JMP IN 

LAW -240 

SZA! I AC 
JMP .-I 

JMS ADCSET 

JMP .+3 

.DSA ONE 

.DSA RETURN 

JMP OFF 

IN LAC (ACFiz 

TCA 

Tad cm 

DAC CUTXR# 

LAC LOIM 

DAC IMPT# 

LAW -740 

DAC M740# 

LAW -4 

DAC MINT# 

DAC STOP 

LAW -240 
SPA! I AC 
JMP .-I 
JMS ADCSET 
JMP .+3 


/FIRST DISK FULL? 
/SECOND DISK FULL? 
/ALL DISKS FULL? 


/ONE sample 


/! 000(10) INPUT DUFFER SIZE 

/ENTRY (LEVEL 7) 

/IGNORE data 
/CHAN 0 NEXT? 


/TRY AGAIN 


/GET RELATIVE DISPLACEMENT 
/TO OUTPUT array 

/RESET INTERMEDIATE ARRAY 

/480 samples PER MIN 

/4 min per outbuff 
/continue data collection 

/waste 250UG 

/RESTART ADC DATA TAKING 


.DSA N 

.DSA PR ADD /ADDRESS OF COHERENT INTEGRATION ROUTINE 

/ ZERO intermediate ARRAY 


law -740 
PAX 
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FSCAT progi'mn (cont.) 


CLl.R 

D/.M 1MaRR-^740.X 

AXS I 

JMP 

/ PLAY IT SAFE-2t:R0 SETUP ACPI?, & MEANS 



LAO 

(JMS 

STORE 


DAG 

PCLR 



DAC 

RCLR 



dac 

I CLP 



I,AC 

(JMS 

MNSTR 


DAC 

MXOLR 



DAC 

MYCLR 


OPI- 

.RI.XIT 

RETURN* 

i 

CLOP 

ISZ 

S 

/TNCRliMENT DAT SLOT 


rsz 

1. 



ISZ 

wate 



ISZ 

PAUS+I 



rsz 

GOO 



ISZ 

CF-RR 

/ALSO PDR TERMlNATinj DUE m 


ISZ 

BDUT 

/IVRITE ON NEXT DISK 


JMP 

S 


CLOO 

.CLOSH 

7 

/ALL DISKS rULI. QUIT 


.EXIT 



PRADO 

a)HINI’+‘)0DW0 


ADERR 

i) 


/ADC FAILURE (CLOCK EXPIRED) 


LAW 

S774 

/ISSUE terminal error 


DAC 

FRCDDE 



I,AC 

aderr 



DAC 

ERA RG 



LAC* 

(?02 



ISA 


/PROTECT THE MONITOR 


JMS 

ERROR 



DBK 




.RLXIT 

ADERR 



• ORPIN 

. I NT M 



LAC 
1 1 S 

FUJFF"I ,X 

fM 


LRSS 

1 V ' 

ID 



TAD 

IT) 1 , X 



DAC 

PDI ,X 



• ENDM 



■.lOHINT 

0 


/ENTRY 


LAC 

riMERR 



SNA 




JMP 

TTT 

/NO ERROR 


I.AC 

COHINT 



DAC 

ERARG 



LAW 

5773 



DAC 

HRCODF 



l-AO* 

(202 



ISA 




rSCAT program (cont.) 


m 


dEGINT 


JMiJ 

DHK 

JMP 

LAC 

rcA 

PAX 

AAC 

PAL 


HRRDR 

mu r 

13PFLO 




. I NT 

0 

. INI' 

5W 

. I NT 

1 zW 

. I NT 

1 70 

,INT 

Z 40 

. I NT 

Jl 0 

. 1 N V 

360 

. I NT 

430 

.INF 

1)00 

. IN r 

‘“>50 

.INT 

6z0 

.INI 

670 

. INT 

740 

. I N I* 

1 01 0 

. INT 

1 060 

. I NT 

1 1 30 

. INT 

1 Z00 

. INT 

1 zS0 

. INT 

1 3z0 

. INT 

1 370 

. INT 

1 440 

. TNT 

1510 

.INT 

! 560 

.INT 

1 630 

.INT 

1700 


SPA 

TAU (-1 

DAC PDl ,X 

AXS I 

JMP B POINT 

/ NOW FORM ACFS 

JMS ACFl z /3I 

/TEST IF MINUTE COMPLETE 
IGZ M740 

JMP QUIT 

/ NEW MIN starts WITH NEXT POINT 


/CONVERT I'D I^S COMP 


6 MSEC 


LAW 

-740 

DAC 

•1740 

LAC 

(JMS 

DAC 

PCLR 

DAC 

RCLR 

DAC 

ICLR 


/RESET minute COUNTER 
STORE 

/SET UP TO CLEAR 
/NEXT ACI* IN 
/UUrPUT BUFFER 
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FSCAT program (cont.) 



JMS 

MNSUB /subtract DC AND CLEAR 


LAC 

QUTXR /POINT TO NEXT ACF 


tad 

N 


DAC 

DUTXR 


lac 

(JMS MNSTR 


DAC 

MXCLR 


DAC 

MYCLR 


ISZ 

MINT 


JMP 

QUIT 

/OUTPUT 

BUFFER 

FULL 


D7.M 

STOP /NO MORE DATA PLEASE 


DZM 

SUBR /IGNORE NEXT BUFFER 


ISZ 

SEO 


i.y\c 

SEQ 


DACi^ 

LAS 

SEQAD 


DAC* 

WDAD 

BOUT 

.REAf.W 

4, CM, 4002., SYNC, 6 


JMP 

QUIT 

/want TO INTERRUPT PROCESSING. FIRST STOP CLOCK 
GLOF 

WATE 

.REALW 

5,4 ,CM,400z,PaUS,0 

QUIT 

LAW 

-1 /FREE BUFFER TO A/D 


DAC 

COMFLG 


.RLXIT 

CDHINT 

FILE 

.SIXBT 

'23MAR' 


EXT .SIXBT '78A' 

/ ACFI2 iNTCGRATES PRESENT DATA WITH rNTERMEDl'ATE ARRAY 
/ AND FORMS DC ESTIMATE. FINALLY PRESENT DATA ARE INCLUDED 


/IN THE 

IMARR 



ACFI 

0 


/ENTRY 


LAC 

BPFLG 



TCA 


/FIND PRESENT DATA 


TAD 

(PDI-I 



DAC* 

(16 

/USE AUTO INC REG 16 FOR X'S 


A AG 

24 

/|7 FOR Y'S 


dac* 

(! 7 



LAC 

IMPT 



DAC* 

(1 4 

/USE 14 & 15 FOR IMARR 


DAC* 

( 15 



TCA 




TAD 

HG IM 



TCA 


/NEC 4 POINTS above IMARR DISCON 1' I NUlTY 


DAC 

GUNTR# 



LAW 

-?.4 

/DO 20 HEIGHTS 


DAC 

CNTHT/^ 



LAC 

OUTXR 



PAX 


/ACPI 2, X REFERENCES OUTPUT ARRAY 

LOnPH 

AAC 

62 

/zb DOUBLES PER HEIGHT 


PAl. 




LAC* 

1 6 

/GET X(T) 
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(''SCAT program (contO 



GSM 


/ABS VALCI^S COMP)iSIGN TO LINK 


DAC 

XYT0 



dac 

XXT0 



DAC 

PX 



LAC 

(CLL 



SZL ! OLL 


/FEST SIGN 


AAC 

2 

/STL-CLL+z 


DAC 

XXT0SN 



XOR 

( z 

/OPPOSITE SIGN FOR KK) 


DAC 

XYT0SM 

/I (K)=YX+(-X)Y 


I.AC 

PX 



MULS 


/CLEAR EAE SIGN?GET X**2 

PX 

0 



PCLR 

JMS 

STORE 

/DR JMS DADD TO SUM 


LAC* 

1 7 

/GET Y(T) 


GSM 


/A3S VALd'S COMP)»STGN TO LINK 


OAC 

YXT0 



DAC 

YYT0 



DAC 

PY 



LAC 

(CLL 



S2L1CLL 


/TEST SIGN 


AAC 

Z 

/STL=C LL-i-2 


DAC 

YYT0SN 



DAC 

YXT0SN 



LAC 

PY 



MULS 


/GET Y**^ 

PY 

0 




JMS 

dadd 

/SUM TO POWER 


AXR 

z 

/POINT TO R( 1 ) 

/ DO R(K) 



LOnPT"XXT0SN 



XXT0SN 

XX 


/STL OR CLL?GET SIGN OF X(T) 


LAC* 

1 4 

/GET Y(T+K) 


MULS 



XXT0 

0 


/FORM PRODUCT 

RCLR 

JMS 

storl: 

,AQR JMS DADD TO SUM 

YYT0SN 

XX 


/STL OR CLL^GET SIGN OF Y(T) 


LAC* 

14 

/GET Y(T+K) 


MULS 


/PRODUCT ( 1 'S COMP) 

YYT0 

0 




JMS 

DADD 

/ADD ru OUTPUT POINT 


axr 


/POINT TO T(K) 

/DO I ( K 

) 



YXT05N 

XX 


/GET SIGN OF Y(T) 


Lac* 

15 

/GET X(T-i-K) 


MULS 


/I COMP PRODUCT 

YXT0 

0 



lOI.R 

JMS 

STOR'; 

/OR JMS DADD TO SUM 

XYT0SN 

XX 


/GFT MINUS SIGN OF X(T) 


LAC-A- 

1 5 

/GET YCT+K) 


MULS 


/| 'S COMP MR 1 Tier 


FSCAT program (cont.) 


XYT0 0 



JMS 

DADD 

/SUM TO I (1C) 


ISZ 

CONTD 

/z POINTS PER LAO 


rsz 

contr 

/TEST FOLDUVER 


JMP 

DOWN 



LAW 

-740 

/FOLD TO BOTTOM OF IMARR 


TAD* 

(15 

/RESET 15 


dag* 

( 1 5 



DAC* 

( 1 4 

/RESET lodlNCE PER ACFI2) 

OOWiJ 

AXS 

2 

/POINT TO NEXT LA'JfTESr aCF DONE 


JMP 

LnOPT 



PXA 




ISZ 

CNTHT 

/TEST LAST HEIGHT 


JMP 

LOOPH 

/DO NEXT .11; I'll IT 

/ ACFS 

DUN El ALWAYS RESET STUR'< I'O DADD 


Lag 

(JMS 

DADD 


DAG 

PCLR 



dac 

RCLR 



DAG 

ICLR 


/ STORE 

PRESENT 

■ DATA 

IN IMARR & SUM(STURR) TO M’-iANS 

/ REVISE IMPT, 

THE IMARR DISCONTINUITY ADORf%S$ 


LAG 

IMPT 



SAD 

LOIM 

/MIN VALUE? 


LAG 

HGIM 

AES, GET .MAX VALUE 


AAC 

-2 

/DECREMENT TO NEW POINT 


DAG 

IMPT 

/SAVE 


tga 




TAD 

LOIM 



TGA 


/# POINTS BELOW DISCONT. 


PAX 




LAG 

(740 



PAL 


/SET MAX M POINTS 


LAW 

-24 



DAG 

CNTHT 

/DO 20 HEIGHTS 


LAG* 

(16 

/1 6 STILL POINTS TO Y^S 


AAG 

-Z4 

/1 7 POINTS TU X'S 


DAG* 

( 17 



LAC 

(MN-I 



DAC* 

( 1 4 

/I4&I5 POINT TO Mil array 


DACw 

(15 

/(DOUBLE WORDS 

LOOP 

LAC* 

17 

/GET X(T) 


DAC 

IMARR 

,X /STORE IN IMARR 

MXCLR 

JMS 

MNSTR 

/OR MNADD TO SU*( TO MN ARRA/ 


LAC* 

1 6 

/GET Y(T) 


DAC 

IMARR 

,X+1 /STJ'JE IN IMARR 

MYCLR 

JMS 

MNSTR 

/STOiiHdJR ADD) TO MM ARRAY 


AXS 

30 

/POUT TO MEXT ilEUHT I -1 


JMP 

DTT 


/ FOLD 

TO NEXT 

HEIGHT 

IN IMARR 


AXR 

-360 



AXR 

-360 

/DECRE»'Ei-IT :n '^Y hU'- 


FSCAT program (cont.) 


DTr 

I H/. 

ci-imr 

/ri-ST LAST HEK^HT 


JMP 

I.llLlP 

/DU NEXT HEIGHT 

/DONE 

WITH PRESENT DAfAlRESHr WNS FR TO MNAHO & RETUR* 


LAC 

(JMS 

MNADD 


oac 

MXCLR 



Dac 

myclr 



JMP* 

ACPI 

/RETURN 

/MNADO CDMVGIM'Jj 

AC TO P'S cusp \ adds to MU ARRAY 

/IN UPi5IDIH)OWN 

FORMAT! I. L. LSWiMSW ) 

MNA'M) 

a 


/RNTRY 


SPA 




I AC 


/CONVERT 7/S COMP 


SHAL 


• /SICN(MSW) TO LINK 


TAD* 

14 

/ADD I.SW'S 


DAC* 

IS 

/HAVE 


SNL 


/TH-JT CARRY XOR SIGN 


JMP 

NOTH 

/PUIl'T ALTER MSW 


SPA! CL A! I AC 

/GET +1 II- PUS 


LAW 

-1 

/GiLT -1 IK IIEO 


TAD/r 

1 4 

/ADD TG -ioW 


DAC* 

1 5 



JMP* 

MNADD 


NOTH 

CLA 


/MU.ST I )CRK’«EUT l4.ilS 


JMP 

.-4 

/anyway 

/ mnsth converts AO ro 

COMP & STORES IN MM ARRAY 

/ (III 

UPJloE-Dtl.'M FURXA D .RHC 14 IS NOT USED. 

MNSli; 

0 


/entry 


SPA 




I AC 


/CONVERT 


DAC* 

IS 

/PUT II' LS.'I 


SPA!CI,A 


/MSW=0 11- PUS 


LAW 

-1 

/Msw=- 1 I F uric 


i)AO* 

IS , 

/PUT IN MSW 


JMP* 

MIloTi? 

/RliTIlRM 

/ siJHPnjrriT to 

SUBTRACT DC(MEANS) F'FlM OUTPUT lATA 

MliSU» 

0 


/ENTRY 


I.AC 

(MN-I 

/I5‘U6 point TD MN ARRAY 


DAC* 

( IS 



DAC* 

( 16 



LAW 

-?4 



DAC 

CNl’HT 

,/DO z0 HEIGHTS 


LAC 

OUTXR 



PAX 


/XR ADDRE'TIES JiJTPiJT ARRA^ 

LOPH 

AAC 

6z 

/z5 DnUBLES f>ER ACF 


PAI. 




JMo 

SQR 

/GET X*X/480 T‘l i)I |/•GIC^:? 


DAC 

TEMPI 

/SAVE DI REGISTE.I 


LAC- 1 




DAC 

TL'MPzy 



JMS 

SOR 

/GET Y*Y/480 IN DI REGS 


FSCAT program (cont.) 


/UUUULF’ ADD(POS DEF SO CAN i'S COMP ADD)8.I'S CO;<P RESULT 


tad 

TEMPI 


DAC 

TEMPI 

/ADD MSW'S 

LAGQ 



CLL 



TAD 

TFMP2 

/ADD LSW'SjCARRY TO LINK 

I.MQ 


/RETURN TO MQ 

LAC 

TEMPI 


SZL 


/TEST CARRY 

I AC 
CMO 

CLLICMA 


/i 'S COMP RESULT 

DAC 

l.ACQ 

TEMPI 

/SAVE FOR LOOP 

DAC 

TEMP2 


LAC 

TEMPI 

/MEAN TO DI REG 

JMS 

DADD 

/ADD (SUBTRACT) TO ZERO LAG( POWER) 

AXR 

2 

/POINT TO R'S 

LAC 

TEMPl 


l.MQ 


/RESTORE -DC TO DI REG 

LAC 

Tempi 


JMS 

DADD 

/ADD TO R(K) 

AXS 

4 

/POINT TO NEXT RiTEST LAST LAG 

JMP 

LOPT 

/DO NEXT lag 

PXA 



ISZ 

CWTHT 

/TEST LAST HEIGHT 

JMP 

LOPH 

/DU NEXT HEIGHT 

JMP^ 

MNSUB 

/RETURN 


/SQR Takes mean, squares it & divides by 4ao 

/TO AVOID OFLO.MEAN IS FIRST DIVIDED BY SRT(480) 
/TO AVOID truncation ERROR, MEAN&SRT( 480) ARE FIRST 
/ADJUSTED TO YEILD AN 18 BIT RESULT THEN SQUARED 
/AND ADJUSTED BACK .( RESULT IS TRUNCATED BY <I ) 


0 


/ENTRY ( <0.1 MSEC) 

LAC* 

l.MQ 

15 

/JET LSWCJPST DE-DOWN FORMAT) 

LAC* 

1 5 

/MSW 

SPA! CLL 


/A0S VaL, CLEAR LINK 

JMS 

C0MP2 

/z^S COMP IF NEG 

DIV 


/DO ROUGH DIVIDE 

25 


/SRT(480) 

SZL 


/TEST OVFLO 

JMP 
CL A 

DCOFLQ 


NORM-1 


/SHIFT TO get 18 BIT RESULT 

LACS 


/GET SC 

Tad 

( 'cZ 


AND 

(37 

/5 BIT SHIFT COUNT 

DAC 

MSHIFT# 

/SAVE 

RCL 


/Times z 

TAD 

(LRS 

0 /FORM IRS ?*M 


/FORM IRS 
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FSCAT program (cont.) 



UAC 

SHFBK /RE-ADJUbT AFTER SQUARE 


LAC 

MSHIFT 


TAD 

(LLS U 


DAC 

.■^6 /LLS M+12 


i.AOA- 

LMO 

16 


LAC-v 

spaicll 

16 /GET MEAN AGAIN 


JMS 

CDMPz /ABS VAL, CLEAR t.INK 


LL.S 

1 4 /LEFT SHIFT P BITS 


DIV 

/DIVIDE BY 


.3 

/SRT( A80)*P'-^’^I 


LACQ 

/GET RESULT 


DAC 

. + z 


MUL 

/SQUARE IT 

SHFBK 

ky 

!.f?3 

0 /SHIFT BACK BY 2M 


JMP-V 

SQR /RETURN 

/t’ilMP? 

performs 

2'S COMP OF Df REG 


0 

/ENTRY 


DAC 

TEMP-^ /ilAVE M.SW 


AcQ 

/GET LSW 


CLUTCA 

LMO 

/COMP LSW 


LAC 

TEMP 


SZL'.CMA 

/CUMP MSW,TEST CARRY 


CLU IAC 

/CARRY IN & CLEAR LINK 


JMP* 

OOMPz /RETURN POS UEF Df 

/DC HF 

LH ICSUES 

ERROR IF DC EXCEDE5 36 BITS 

DOUI'f.lJ 

I.AC'.^ 

(Z0Z 


ISA 

/PROTECT MONITCR 


LAC 

COHiNT 


DAC 

RRARG 


Law 

S772 


DAC 

ercdde 


.IMS 

DB< 

ERROR 


JMP 

QUIT /EXIT PGM 

/Ai'lUHi: 

STORES THE DI REG IN I'HE OUTPUT RUFF 

Sl'iJFH 

0 

/ENTRY 


DAC 

LACQ 

ACFI2,X /STORE MSW 


UAC 

ACFIz,X+l /STORE LSW 


JMP^ 

STORE /RETURN 

/UAUfJ 

adds ( I 

COMP) DI REG TO OUTPUT BUFF 

,)AIjU 

0 

DLL 

/ENTRY 


TAD 

ACFIz,X /ADD MSW^SsCARRY TO LINK 


DAC 

ACFI2,X /STORE 


LACQ 

SZU.CLL 


/TEST I'AliPY 


FSCAT program (cont.) 



I AC 


/CRRY AROUNDiCARRY OUT TO LINK 


TAD 

ACFlz,X+l /ADD LSW/SiCARRY OUT TO 


SNU 


/TEST carry 


JNP 

OUT 



isz 

ACFI 2.,X 

/CARRY INfSKiP ON CARRY AROUND 


SKP 




lAC 


/carry around 

OUT 

DAC 

ACFU.X*I /SAVE LSW 


JMP* 

DADD 

/return 

A DC SET 

0 

JMS* 

.DA 

/ENTRY 


JMP 

PAST 


WC 

. DSA 

WC 

/ADDRESS OF WORD COUNT 

SUBR 

.DSA 

SUBR 

/ADDRESS OF 'SUBR' ARG. 

PAST 

LAC 

( 404000 

/RAISE TO PRIORITY LEVEL 4 


ISA 

LAC-A 

( 1 55 

/SETUP ADC lONCE ONLY 

REAL 

DAC 



COMFLG 

JMS* 

.-1 


DIFPC 

703701 



BPFLG 

ADS VC 
DBK 



SVAC 

LAC* 

( 151 


ERADD 

DAC 

REAL 



LAC 

(JMP 

INI 


DAC 

PAST 


INI 

LAW 

-1 



DAC 

COMFLG 



TAD 

( BUFFI 



DAC 

BUFI# 



LAW 

-1 



TAD 

(BUFF2. 



DAC 

BUF'^# 



TCA 




TAD 

BUFI 



DAC 

DIFFC 



DAC 

BPFLG 



LAC* 

TCA 

WC 



DAC 

WC 



LAC* 

SUBR 



DAC 

SUBR 



LAC* 

(Z0Z 



ISA 




JMS 

ADIN 



DBK 




JMP* 

adcset 

/RETURN 

ADIN 

0 


/ENTRYY 


LAC CL0C<+3 /RESET CLOCK 

DAC.v (7 

LAC WC 


/SETUP DCH 


FSCAT program (cont.) 


ADSVC 


ADG 


EXIT 


SLOWP 


ERROR 


DACv 

(?.6 


LAC 

BPFLG 


TAD 

BUFz 


DAC* 

703704 

(z7 


703744 


/clear flags 

703724 


/ENABLE TRANSFERS 

JMP* 

AOIN 

/RETURN 

0 


/LEVEL 0 ENTRY 

DAC 

SVAC 

/SAVE AC 

703704 
LAC 
I AC 
5ZA 

M740 

/CLEAR OFLO 

JMP 

ADG 


LAC 

MINT 


I AC 



S7.A 



SKP 



DZM 

STOP 

/DISABLE A/D DURING BUFFEROUT 

703Vzl 


/TIMING ERROR?? 

JMP 

.+3 

/NO 

ISZ 

Tim ERR 

/INFORM USER 

703744 


/CL AAR FLAG 

rsz 

CDMFLG 

/set busy flg zero 

JMP 

SLDWP 

/IT ALREADY WAS— ERROR 

LAC 

SUBR 


SZA 



JMS* 

REAL 

/PRIME TO RUN SUBK AFTER EXIT 

LAC 

DIFFC 

/SWITCH BUFFERS 

XOR 

BPFI.G 


DAC 

BPFLG 


LAC 

STOP 


5ZA 

JMS 

AD IN 


l-AC 

( 404000 

/REQUEST monitor AFTER EXIT 

ISA 



LAC 

SVAC 


DBR 

JMP* 

ADSVC 


LAC 

ADSVC 


DAC 

ERARG 


LAW 

5777 


DAC 

ERCODE 


IMS 

ERROR 


JMP 

EXIT 


0 


/entry TD error PRINTOUT 

LAC* 

( 166 

-T^ONITOR ERROR SUBROUTINE 

DAC 

ERADD 


LAW 

5770 


JMS* 

ERADD 



ERG ODE 


FSCAT program (cent.) 


f^nARG 

mni 

DZM 

STOP 

/ARGUMEHT 

CERR 

.CLOSE 

5 



JMP* 

ERROR 

/RETURN 

BUF(“ 1 

. BLOCK 

17B0 


PDI 

. BLOCK 

50 


BUFF 2 

. BLOCK 

1750 


PD2 

.BLOCK 

50 


IMARR 

. BLOCK 

740 


MN 

. BLOCK 

1 20 


LOIM 

.DSA 

IMARR-i 


HGIM 

.pSA 

.EMD 

IMARR+7 37 
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APPBNDIX IV POSl’-PHOCRSSlNG PROGRAMS 
Th<i programs listed below are for post-procossing of the correlation 
functions which were obtained in real time and stored on disk or OCCtnpe, 
the source files for the programs POWW, VELh, AREAD, and DINFLT are loaded 
onto disk 3 of the POP- 15 from the program DBGtape, the editor program of 
the system software is used to modify the file pointer in the ARRAO program 
so as to agree with the file and extension of the first file to be procossedt 
the statements are accessed with the edit command "L jSIX". After closing 
the modified ARBAD source file the MACRO compiler is used to produce binary 
files for these two programs, the POOTRAN compiler is used for POWW and 
VELL, Bach time a new data set is processed AHEAD must be edited and com- 
piled with the now file name. The source of the data, either disk or DRCtape 
is assigned to ,DAt slot 5, the teletype to slot 2, the paper punch to slot 
3 and disk 3 to slot -4, The programs POWW, AHEAD, and DINFbT are then loaded 
and executed to produce a paper tape of power data, Simiiafly vBUj, AkHAD, 
and DINFLT are used to produce a velocity paper tape. 

AHEAD program 



.GLOBL 

set, GET,. DA, DINFLT 


. lODEV 

5 

SET 

0 

/ENTRY DIMF.N. 


. IN IT 

5,0, SET 


.SEEK 

5 , F I LE 


ISZ 

EXT 


JMP* 

SET 

FILE 

.SIXBT 

'04APR' 

EXT 

.SIXBT 

'V8A' 

GET 

0 

/ENTRY 


JMS* 

.DA 


JMP 

. fz 

ARH 

0 



LAC^ 

ARH 


dac 

R>2 


dag 

arg 

R 

♦ READ 

5,4,ARR,4002 


.WAIT 

5 


J MS ★ 

DINFLT 


JMP 

,*z 

ARG 

0 



JMP^ 

.END 

GEr 


DINPi;r progrom 


iJlNFLr 


\RO 


POS 


PAST 


SIGN 


.TITLE 

DINFLT 


.ULOBL 

.da.dini-lt 

0 


/I'NTRY 

JMS* 

.DA 

/GET ARG 

JMP 

. '^ Z 


XX 



LAC 

(ARC 

/ARRAY ADDRESS 

TCA 



TAD 

ARO 

/RELATTVE ADDRESS 

PAX 


/USE XR 

TAD 

(7640 


Pal 


/USB LR 

LAC 

ARG.X 

/GET MSW 

SMAICLL 


/CHECK SIGN 

JMP 

PUS 


LAC 

( 400000 


DAC 

SIGN 

/SET sign NEC. 

LAC 

ARe,X«H 

/GET LSW 

CMA 


/ABS VALv OVFLD TO LINK 

LMQ 



LAC 

ARO\X 

/get MSW 

CMA 


/ABS YALl TEST LINK 

JMP 

PAST 


DZM 

SIGN 

ASET SIGN POS. 

Lac 

ARG.X+I 

A5ET LSW 

LMQ 



LAC 

ARG.X 

/GET MSW 

NURM-I 


/LEFT SHIFT BOTH 

XOR 

SIGN 

/attach SIGN 

DAC 

ARG,X+ 1 

/PUT IN SECOND WORD 

LAOO 



AND 

(777000 

/CLEAR EXP BITS 

DSC 


/ATTACH EXP 

SNA 


/IGNORE ZERO 

JMP 

. +3 


XOR 

(77 


I AC 



DAC 

ARG.X 

/PUT IN FIRST WORD 

AXS 

z 

/NEXT PAIR 

JMP 

LOOP 


JMP-A 

1 A 

DINPLT 

/RF.TURN 

.END 
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POWW program 


DIMENSION IN(4002) 

COMMON B(25,20,4),P(20) 
equivalence (IN(I ).B(I ,1.1 )) 

PMAX=0. 

NR=0 

DO 700 IFIL=I , I 0 
CALL SET 
DO 700 IREG=| ,3 
CALL GET (IN) 

DO 100 I M= I , 4 
NR=NR+| 

DO 200 [H=|,z0 

IF(B(I , IH, IM).GT. 1. )P( IH)=ALOGIO(B( 1 , TH, IW)) 
IF(P( IH).GT.PMAX)PMAX==P(IH) 

200 CONTINUE 

REFLOAT (NR) 

WRITE (3, 300 )R, (P( IH) , IH=| ,20) 

300 FORMAT(F5.0,20F5.2) 

100 CONTINUE 

700 CONTINUE 

WRITE (2, 340) PMAX 

340 FORMAT! 1 H0,6HPMAX :=,FI0.3) 

STOP 

END 


VELL program 


dimension IN(4002) 

COMMON B(25,20,4),AMP(20,I 3),RF(z 0. 13).AIM(20, 13),V(?0,4) 
COMMUN/PATA A/PT( 20, 1 z8) , ALT( z0) 

EQUIVALENCE ( IN(| ) ,B( 1 , 1 , I ) ) 

VMAX=('i. 

NR-0 

DO I 0 1 = 1 , z0 
ALT(I)=60. + I .5'VFLQAT(t) 

DO i 0 J.= I 2 I , I 20 
10 VPT(I,J)=0. 

DO 700 IFIL= 1.10 
CALL SET 
DO 700 IREC=I ,3 
CALI. GET(IN) 

DO 100 IM=1,4 
NR=NR+1 

DO 200 IH=1,Z0 
L= 1 

RE( IH.L)=B( 1 , IH, IM) 

DO 50 J=2.z5,z 
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Vlit.L program (cont.) 


L=L+ 1 

RE(IH,L)sH(J,IH, IM) 

50 AIM(IH,L)-B(J+I , IH, IM) 

^00 CtlNTINUK 

DO 2.07. I H=l , 7.0 
DC) Z07. IL-I ,4 
7.07, V(TH, rL)-“Q'. 

DO 290 IH==| ,20 
DU 270 IL-2,4 

AMPCIH, IL):=SQ«T(RECTH, I L) **2+AIM( CH, rL)-.»(^2) 
fF<AMP(IH,IU.LT.(.|0*RE(TH,l )))GU TQ ?62 
IFCREC IH, lU . GQ.0. )CO TO 260 

V ( r H , r L ) ^4 . 67MTAN2 ( A TM ( I H, I L ) , RE ( TH , f L ) ) /FL0AT( F L- I ) 
GO TO ?,70 

260 V(IH, IL)=’7.335/FLaAT(TL-l ) 

GU TO ?,7.0 

262 V(lH,rL)"0. 

I.MAX = 1L“I 
GO TO 7.75 
270 CONTINUE 

LMAX=‘i 
275 VAI-0. 

VA2«0. 

DO 280 I L= 2, UMAX 
VA 1~VA l+AMP(IH, IL)*V( IH, ID 
280 VA^-VAz+AMPC IH, ID 

VPT( IMyNR)=VAl /VA2 
ABV»:ABS (VPT(IH,NR)) 

IF(ABV.GT.VMAX)VMAX=ABV 

IF(AbV.GT,9.99) VPT (I H, NR ) =S ICN (9. 99,VPT( IH,NR) ) 

2.90 CONTINUE 

100 CONTINUE 

700 CONTINUE 

DO 300 IH=I,20 
WRITE(3,333)ALTf IH) 

333 FORMAT ( IH ,F5. I ) 

ND=1 
N2-,= 16 

DO 340 NNt=l,8 

!'VRITE(3, 344) (YPT(TH,NR ) ,NR=NI .Nr) 

344 FORMATCIH ,16F5.7.) 

N! -Nl -H 6 
Nr-Nr+i 6 
340 CONTINUE 

300 CONTINUE 

WRITE(7., 440.WMAX 

440 FORMAT ( I H0 , 6HVMAX “ . F 1 0. 3 ) 

500 STOP 

END 
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APPRNDIX V PLOTTING PROGIUMS 

The following programs are run on a Hewlett- Packard 9830A computer with 
16 K words of memory. The computer is equipped with a thermal printer, paper 
tape reader and a plotter. A cassette drive is an Integral part of the 
machine. 

The first two programs listed below arc used to read in a p.aper tape 
and store the data in a more compact form on cassette, liach data cassette 
is marked with files numbei-ed 0 to 11 of length 2770 words. The storage 
programs: 1) ask for file identification data, 2) read the paper tape and 

convert to integers, and 3) store the data in a user-specified file on the 
data tape. As the file identification information is input the printer 
produces a hard copy which is stored with the cassette in a filing system. 

The final three programs plot the data stored on cassette. The user 
inputs a file number and the computer reads the data from tape and indicates 
what file is about to be plotted. The programs will supply recommended plotter 
control parameters but the user can override and choose different values. 
Generally a plot is made using the recommended values and if unusual signal 
characteristics occur then a second plot is made with new plot parameters. 

The programs can loop back and produce any number of plots witliout re-reading 
the data from cassette. 

1, Velocity data storage program 


111 . n iiu l;j8.n Lfl •nj].H3.M3»H'3»S2(L2»fl6ifl?>S4)H2)El>E2s(-3.<Ea.E5 
:-v Mil Vii, MORE VER'inn 11»06»06»23»78 IN»PfiPEP TRPE.OUTBCftSSETTE 
.'1 MM Mnltl DU enSSETTt; IM IMT. PPPRY D. RQIU-U01‘) HT. B!l! 
ill ViitI hU i.2H:PWSCai.> 1.^83 
4ii I t.f I 1 TO 128 
Ml Hi I IM 
I..II III.: I 1 
ru III 1 M 'iq 

DO I'll 
■•VI )l[ 

S‘r> pr-itir 

IDO ^■‘RINT "DflTC?- 
lie JMPUT 1 $ 

J20 PF'IMT l.r 

I'iO u[LPn;u—n:i=‘' velocities 

Mfi RRiHT "STfiKT TIME'’ INPUT HOURSMIIHUTES. “ 
l‘;0 lill’UT h4,M3 

1, 'U fRTlir 'DflYLIGHT SlTVIHCS TIME?" 
toil I'l lNT INPUT 6 FOR CSTjS FOR CDT" 

I'^e IN' Ui 29 

r»l It E9-6 

I'M i UK ■•■STHFT«"iH l^lOO+MS! "CST" 

.•JII ll'im 'HO. OF RECORDS'?" 

2, ’0 nil 'LIT N3 
.'.111 11:1111 H3 


2-I0 PPINT 'HS-EflSE ALT, El-STflRT SETTING^ INPUT El" 

2Al IMPIJT El 

242 H.-'==CE1“2:.»1.^ 

251 PRINT ''H2-"JH2 
25S PRINT "VELOCITY" 

2f’0 rOR 12=1 TO 20 
280 ENTER a,290ffl)H8 
280 FORMAT 1X»F5.1.5X 
300 FOR 11=1 TO 8 

310 ENTER <.1.320»fl)<F0RJ=lT016iWCI2»<Il-l)*16+J]) 
320 FORMFIT 1X.I3F5.2.3X 
330 NEXT II 
3‘IG NEXT 12 

S50 REM INPUT TONE » CONVERT TO IHT. 

Sm FOR n«l TO H3 

3?0 FOR 12=1 TO 20 

380 liU?,! in«wn2>mrioo 

390 NEXT 12 

400 NEXT II 

405 REM ‘..TORE 

410 PRINT "PLACE DATA TAPE IN TAPE DRIVE" 

420 PPIlir "FILE HO. WHERE DATA IS TO 0E STORED?" 

430 INPUT FI 
440 PRINT "F1LE"IF'1 
450 STORE DATA FI 
480 STOP 
470 Etin 


2. Power data storage program 


IK ri.in nn .i'lf l28T(L<C4tn,H3»M3»N3»S2»L2fA6ffl7»24»H2,El.E2fE3fE4>E5 

2H KEM l■■llNL‘R STORE J 0: 42: 06« 23« 78 IN=PAPER TAPE» OUT=CRSSETTE 

21 Pel' Uirm OH cassette in INT. array D. ROW1=LOH NT. BIN 

3U HIM HII 128]. WSC21> 1283 

40 FOR I -1 TO 128 

50 nri;i-i 

80 tlEXl I 

70 mCIOT^kT 

80 Ari0]-”32 

90 AC 1 2 ]wT2 

95 PRINT 

100 PPIin "DATE?" 

110 INPUT L# 
i2G PRINT LT- 

130 L:rCLEIKLt.’> + l ]=“ POWER CLOG PLOT)" 

MO PRINT "START TIME';> INPUT HOURS. MINUTES. " 

150 INPUT H9.M3 

170 PRINT "DAYLIGHT SAVINGS TIME?" 

180 PRIMT "INPUT 6 FOR CST.5 FOR CDT" 

190 INPUT Z9 

191 H3=H9+29-6 

193 PRINT "START="?H3^^100+N3»'‘C3T" 

210 PRINT "NO. OF RECORDS?" 

220 INPUT N3 
230 PRINT N3 

240 PRINT "H 2 =eASE ALT. E1=START SETTING. INPUT El" 

241 INPUT £1 

242 H2=v:E1“2)»1.5 
251 PRINT "H2="IH2 
255 FOR 1=1 TO N3 

25S ENTER < 1 . 261 . A)WC 21 . I ]. <F0RJ=1 TO20. WC J. I 3) 

261 FORMAT F5. 8. 20F5. 2. 4X 

264 NEXT I 

267 REN INPUT DONE 

270 REM MIN..MAX.AVE 

273 A5=0 

276 S2=15 

279 L2=0 

282 FOR J=1 TO N3 


235 FOR TO 20 

233 IF WCbJ3>6.75 THEM 3C.3 

291 IF JWl THEN 300 

294 I'lC I > J I«7 

297 flOTO 303 

300 Hni.j>wniJ-n 

303 IF WCIiJ] <« L2 THEN 309 

303 L2»Wri!J1 

309 IF WCI»J] S2 THEN 315 
312 S2sNCI,J] 

315 fl5=>fl5+WC Ii Jl 
318 NENT I 
321 NEXT J 
324 flt>Ft5.‘'C20*H3) 

327 PRINT "POMER" 

330 WRITE aSi 339)32, L2»fl6 

333 PRINT "HfiX IS"M0)KL2-S2) 1"DB ABOVE MIN.” 

33CI PRINT "AVE !S“} ie*‘A6-S2)J "DB ABOVE MIN. " 

339 1" OPMHT " M I " » F5 . 2 , 3X , MAX» " , F6. 2 , 3X , " AVE = " , F6 . 2 
350 PGM INPUT DONE, CONVERT TO INT. 

3r.O l-OF' MbI to N3 
370 FOF- I2«l TO 20 
380 nn2, II ]=wi 12, m^ioo 
390 tU*::T 12 
400 NFMT II 

410 PRINT "PLACE DATA TAPE IN TAPE DRIVE" 

430 PRINT "FILE NO. WHERE DATA IS TO BE STORED?" 

430 iiirUT FI 
440 PRINT ■FILE"fFl 
4Uu ^>10PE DATA FI 
Am STOP 
470 CNIi 


3. Velocity data plot program 


Ji,i I i 4 | Iiu 21, 12S'J,L*r.40],H3,M3,N3,S2,L2,A6,A7.24,H2,El,E2,E3,E4,E5 
;*i. i Ut vr;i PLOT VERSION 21 : 40: 10t 31 1 78 
14 I, 

4M liIM HU 128I,Tti;3 3 

50 PPIHT "PLACE DATA TAPE IN TAPE DRIVE" 

CO PF'llir "FILE NO. WHERE DATA IS STORED?" 

70 KIF'IJI PI 
m PRINT "FILE"! FI 
90 LDAI'i DATA FI 
lOU PR i N't L* 

120 T;t‘-"i:ST" 

150 PRINT "LOCRL START T IME" ? M3+H3*-100! T* 

160 PRINT "HO. OF RECORDS=" ! N3 
170 PRINT "H'2=EflSE ALT.="!H2 

190 PRINT "RECOMMENDED PLOT PARAMETERS? 1*YES,0=NO" 

200 INPUT R9 

•210' R9+1 OF 1100,1020 
220 PR 'INT "SCALE FOR 1 M/'S-"JH 
230 PRINT "PEN UP LIMIT M-'S="L3 

240 Nl=-H2 

24 1 T6b5- < MS /S- 1 NT C M3/5 ) ) 

242 T7-- 1S*<. MS' 1 5- 1 NT < M3/ 15)) 

243 IF Tb'>D THEN £45 

244 T6b5 

245 IF T?;0 THEN 247 

246 V7--=15 

247 REM T6.T7 DONE 
250 REM AXIS 

260 SCALE -80,700,-15,115 
270 XR'/.IS 0,-l,0, <6-T6)*5 
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375 0»25»<6-T6)#5»<6-T6)*5+<.INT<N3/'9>>*25 

Z?C IPLOT lOfOil 

277 IfiBEl 1.2» 1.7f0f8/lOTI 

2'7n VflIUS 0*5,Cfl.OS 

390 i05»-if 8f <s-T6)#5 

2'35 liflNlS 105»25t <6-T6)*5» <6-T6)*S*<INT<N3/5>>*25 

300 PEM HEIGHT LABELS 

310 PLOT 0i0»i 

320 FOR I»1 TO 20 STEP 2 

330 IPLOT OflOfl 

340 LABEL <*» 1 . 5» 1 . 7t 0» 8/1 1 > 

350 CPLOT -Vi - 0.3 
360 Hl«HH-3 

370 LABEL <380i 1 . 5i 1 . 7i 0» 8/1 DHl 

380 FORMAT F6.1 

330 HEMT I 

400 PEM TITLE 

410 IPLOT IOO 1 -II 2 .I 

420 LABEL <»i 1 . S. 1 . 7i 0i 8/1 1 )L# 

430 REM JJAXIS LABEL 

440 FOP I5?»16-T7 TO M3 STEP 15 

450 M4SM3+I5-1 

470 PLOT 5*15. Ill 

500 T0« 1 OO*H3+M4+40* I HT < M4/80) 

510 LABEL C*. 1.5. 1.7*e> 

520 CPLOT -0.3. -0.5 

530 LABEL ^540. 1 . 2. 1 . 7. 270. 8/1 1 )TB 

540 FAPMAT F5.0 

550 HEHT 15 

560 REM AH IS DONE 

570 FOR 11“1 TO 20 

580 OFFSET 0.5*11 

590 XAXIS 0. 1.650-<128-N3>*5.0 

600 REM C2— 1 FOR LAST ZERO 

610 REM C2-1 FOR LAST HOT ZERO 

620 C2=-l 

630 FOR 12-1 TO N3 
640 X«I2*5 

650 Y-DCIl. I2]«H/100 , , 

• 660 IF Y=0 OR ABS<Y) >- L3*H THEN 710 
670 REM NOT ZERO 
680 Cl =2 
690 GOTO 720 
700 REM ZERO 
710 Cl-1 
720 C3=C1*C2 
730 C2=<-l>tCl 
740 PLOT X.Y.Cc 
750 NEXT 12 
760 FEH 
770 NEXT 11 
7S0 PEM SCALE 
790 OFFSET 655.55 
800 FOR 1=1 TO 109 
810 B4=I*H 

820 IF B4 >■= 5 THEN 840 
830 NEXT I 

840 YRXIS 0,.B4.-1*B4.B4 
850 LABEL <*.1.2.1.7.0.8/11) 

860 FOR 1=1 TO 3 
870 Ul=B4*<I-2> 

880 PLOT 15.U1.1 

890 CPLOT 0.3. -0.3 

900 LABEL <910. 1.2. 1.7.0»8/I1)U1/H 

910 FORMAT F3.0 

920 NEXT 1 

930 CPLOT -3.2 

940 LABEL f*» 1 . 2. I . 7. 0. 8/1 1 ) "SCALE" 

950 PRINT "HEW PLOT WITH THIS DATA? PREPARE PLOTTER 
960 PRINT "1=YES" 

970 INPUT R 

980 IF R=1 THEN 190 
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990 PRINT 
1000 STOP 

1010 REH PLOT FRRfl. ROUTINES 

1020 IF H2>36 THEN lOEO 

1030 H«5 

1040 L3>1 

1050 GOTO 1080 

1060 H»1 

10^0 L3»7.5 

1080 RETURN 

1090 STOP 

1091 REti VERT. TICKSsS PLOTTER UNITS. H*NO. OF PLOTTER UNITS FOR 1 M/S 
1100 PRINT '‘SCPLE FOR 1 M/S?" 

1110 INPUT H 

1111 REM PBS. VPLUE OF VEL. < L3 TO BE PLOTTED. 

1120 FRUIT "PEN UP LIMIT M/S?" 

1130 INPUT L3 
1140 RETURN 
1150 STOP 
1160 END 


*1. Power versus time at a fixed altitude plot program 

lU LHM im 21* 120;itL»C40]iH3»M3»N3»S2.L2tfi6>Rr»24»H2>EljE2»E3»E4»E5 
?0 PP'I PMIIER PLOT H0FI2. 14i00:0?i 15t?8 
30 Btw 

4U DIM T.'n283»TfC3] 

50 PFIIIT "PLfiCE HflTfl TAPE IN TAPE DRIVE" 
t’M PFINT "FILE HCi. WHERE DATA IS STORED?" 
ru lirUT FI 
80 PPINT "FILE"1F1 

9U unrin data fi 

Ifii't PPINT U 

"r 

15U PPINT "LPCPL START TIHE"jH3+H3irl0OJT4 
16u PPINT "HO. OF REG0 RDSb"sN 3 
IFl llPin: aS) 164)S2»L2»B6 

ItV I'PIIIT "MAN IS"; 10?<L2-S2>;"DB ABOVE MIN." 

U'3 PPINT "AVE IS": 10*tA6-S2); "DB ABOVE MIN." 

164 FuPMiTr "MIh="»F5.2,3K> "NAX=" » F6. 2» 3N» "AVE="»F6.2 
lUS S' ALE 'SOj ? 00i -15' 1 15 
1 /n PPim ‘H2SEASE ALT.-":H 2 

l'n:i PPINT "RECOMMENDED PLOT PARAMETERS? 1=VES>0-N0" 

2A0 INPUT P9 

.•■lO i .r.iiF pq+l 0F-1?2D»1830 

212 PPINT "BAi.Ei="B) "HEIGHT="H.< "LIHIT="L3 

215 PEN THRESHOLD IS T5 

2.,‘U T5''U5 10+S2 

222 UFKSET 0.0 

k* k* fm*- PLOT 5. 111.1 

224 LhUEI, ':225fl.2.1.7.e.e,'inS2..L2.e>L3 

225 FOPMAT " M 1 1"= " . F5 . 2 , 3N . " M AH= " . F6 . 2 . 3X . " BASE= " , FS. 2 . 3X . "L I M I T=” , F3 . 0 » X. " DB " 

226 FEN 

240 H1=H2 

241 TC»5^..M3/S-IHT<M3/5)) 

242 T r- 1 U- • M3 1 S- 1 NT C M3/ 155) 

243 II- U:. 0 THEN 245 

245 IF T7 0 THEN 247 
5-lr roHE 

m\ H.;is 

2r0 MHKIS CI»-1 J0J <6-16 j*' 

::n:;rr* <6-T6 ^^5, ^-T6)*5+aNT<N3/5))*25 

I PLOT 

2r? LABEL 1.2? 

280 ta::i? o.s.o.ios 

290 MAX IS 1 05 . - 1 ! 0 , < c-T6 5 *5 


295 105»£5. C6“T6>*5+<1NT<H3/'3>)«25 

300 PEtl HEIGHT LABELS 

310 PLOT 0<0fl 

320 FOP I«1 TO 20 STEP 2 

330 I PLOT Oil0»l 

340 LfiDEl. 1 .5» i.?» 0 .e.'lO 

350 CPLOT -7, -0.3 

seo HtBHl+3 

370 LABEL i 380» 1 . 5» 1 . 7» 0» 8/1 1 )Hi 

•330 POPNflT F6.1 

390 RLLT I 

400 PLH TITLE 

410 ir-lOT 100,-1 12il 

4.:0 label f«,1.5, 1.7(0,8 llJLf 

4 in PLH I'fl'IlS LABEL 

440 fOP i5«lg-T7 TO 113 STEP 15 

4‘T0 

■L.’u MUT 5M5,1,1 

r.011 Tf.= 1 00 * H3+P4+4G* 1 1 IT ( M4 -80 ) 

510 LABEL 1.5, 1.7,0) 

520 CPLOT -0.3, -0.5 

530 LABEL CS40, 1 . 2, 1 . 270, 8/1 1 )T6 

540 FORMAT F5.0 

550 me::t is 

5G0 REM AMIS DONE 

920 REM GIVEN HEIGHT VS TIME 

930 REM IHIT. T ARRAY 

940 FOR I>1 TO 128 

950 TCKIbO 

980 NE,‘!T K 

l':40 FOP H9sl TO 20 

1350 OFPCST G?£i(H9 

1355 C2P-1 

1500 FOP IlBl TO H3 

1370 :;sii#5 

13G0 I2«H9-1 

1390 C'OSUB 1900 

1400 V=<2n'8)*H 

1410 G3SUE 1000 

1420 HEP.T !1 

1430 PEN 

1440 HEKT H9 

1470 PEN 

1510 REM LABEL 10 DB SCALE 
1520 OFFSET 0,0 
1530 YRXIS 6S5,H, 105-H, 105 
1540 LABEL <*,1.2,1.7,0,8/11) 

1550 CPLOT -2.£,0.3 

1560 LABEL <*, 1 . 2, 1 - 7 , 0, B-'l 1 ) " 10 DB" 

1562 PRINT "HEM PLOT WITH THIS DATA? PREPARE PLOTTER" 

1564 PRINT "1=YES" 

1566 INPUT R 
1568 IF P--1 THEN 190 
1570 PRINT 
1572 STOP 

1580 i-EM t 2--1 FOR LhST HIDDEN 

IT/'O PPM 1-2=1 FOR LAST NOT HIDDEN 

IvOO IF •TCm-5)>V THEN 1660 

1610 PEM HOT HIDDEN 

1620 TCI1]»V 

16';0 Cl =2 

1640 iTOTO 1680 

1650 PEM HIDDEN 

1660 TCn ]=TCin-5 

1670 01=1 

16S0 C 3=01-^02 

1690 C2=t-1 ,)tCl 

1700 PLOT M,Y,C3 

1710 RETURN 

1715 FEM height FACTOR H DETERMINES THE HU. Of VERTICAL PLOT UNITS 

1716 PEM FUR 10 DB. NOTE THAT VERT. TlCIiS APE 5 PLOT UNITS. 

1720 PPim "INPUT HEIGHT FACTOR" 

1730 INPUT H 

1.-40 PEM BASE VALUE B IS SUBTRACTED FROM ALL DATA BEFORE PLOTTING. 
1.’'41 REM B SHOLLD BE NEAR THE HIM. AND SLIGHTLY HIGHER. 
l.'5U PRINT " INPUT BA..E VALUE" 
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i^:o ^npur B 

l?TQ r:;M LIMIT VFilLUE IS THE LIMIT III DB FiEOVE THE MIIL hBOVE 
irri RLM UHICH THE PLOT HILL BE CLIPPED IF fi i Mill, SPII L nccupg* 
irci? PPIHT "IHFUT LIMIT VRLUE III DB" 

1T90 INPUT L3 
urn L3*IHTa3’> 

IPlCi PEIUPM 
1C20 :.TOF 
1U30 H«IO 
lUlO 45 

inro L 

im pv^um 

IP?0 ’.TOP 

lUC'D rCM LIMITING 

l?10 if 1 1-113 THEN lU4u 

lULO IF TS OR DCIOrl, Il + i ] 100 T5 THEN 1S4U 

is::o ;S:^T5 
1^^4111 rLTULN 
\m END 


5» Power vor^^us height at a fixed time plot program 


IN 

,30 

so 

40 
50 
SO 
70 
Co 
SO 
lOG 
120 
150 
1 BO 
161 
Uv! 

163 

164 
168 
170 
130 
200 
210 
212 
21S 
220 
222 

22^ 

225 

2s6 

240 

241 

242 
24*3 

244 

245 

246 

247 
250 
270 


1 23 J» U*MO 3. H3» M3» N3» S2» L2. R6» fl7i 24» Hgf Ei I E2i E3» E4 1 E5 
PLOT VEPT, 301231 1H01I78 


I'uM iiir^i 
RE'I PII13ER 

ncG 

HIM T')C128 3)T«3 3 

PPIIIT ''PinCE nflTfl TAPE IN TAPE DRIVE" 
PPIHT "FILE HO. WHERE DATA IS STORED?" 
IHPIJT FI 
PPIHT ■■FILE"IF1 
LOAD DATA FI 
PRINT L:T 
Tj=‘‘C'IT" 

PRIlir "LOCPL START TIME" I M3+H3»100( T* 

PRINT "NO. OF REC0RDSt="ni3 

WRITE (15) 1S43S2,L2»A6 

PRINT "MAX IS"no^<:L 2 -S 2 )i "DB ABOVE MIH. " 

PRINT "AVE IS"! 10ir<A6-S2)r‘DB ABOVE MIN." 

FORMAT "MIN“")F5.2)3X)"MAX="»F6.2)3X» "AVE 


'F6.2 


1«VES)0=NO" 


SCALE -80)700.-15,115 
PRINT “H2«*E’ASE ALT.«"!H2 
PRINT "RECOMMENDED PLOT PARAMETERS'? 

INPUT R9 
GOSUB R'9+1 OF 1720,1830 
PRINT "BA'3E="B, "HEICHT="H, "LIM1T'="L3 
REM THRESHOLD IS T5 
T5-L3 '10+S2 
OFFSET 0,0 
PLOT 5,111,1 

LABEL '225, 1.2, 1.7,0, 8/ll)S2»L2,B,L3 

FORMAT ”HIN»",P5.2,3X, "MAN«" , F6. 2, 3X, "BASE=" , F5. 2, 3R, "LIMIT"" , F3. 0, N, "DB" 

PEN 

H1=H2 

TE«5*' N3 5-INTcM3.'53> 

T?‘*. 1 ( M3 1 5- 1 HT < M3/ 1 5 > > 

IF TE 0 THEN 245 
Tb-S 

IP I.' A THEM 247 
T7-n'. 


PLN 

FFN 


fr,T7 rONE 

a::is 


5 ;;) #25 


, “.'e 

23'*; 

200 


::b;:i.s o>-i,0i<6-Tb;i#5 
::m ;l* Q,25* CtW6*^*5; ^6-Tb)*5+aHT« lie 
lOsCbi 

LhFLL 1.7j0j8/imi^ 

I Hi T:’ D » 5> G » 1 05 

NHL* I S 1 05 ! - 1 1 0 ? ( d—T 6)^5 

iq5jS5» <6~tb>* 5» <6-T6)^5+aiiT<N3p-5))^^25 
FHM HLIGHT LBBELS 








Y\ 


# 


r 
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JIU F’LiiT U.U>1 

FOP I I TO 20 STEP 2 
..JJU I plot O»l0fl 
m LUPCi. i.5i 1,7. o» evil J 
350 CPLOT -7»“0.3 
300 Hl»Hl+3 

370 LABEL <330» 1 . 5» 1 . 7» 3i 8/1 1 >H1 
300 FORMAT F6. 1 

390 he::t i 

40U PEM TITLE 

410 IPLOT 100.-112.1 

4^:0 LABEL <». 1.5. 1.7.0. 8/ll)L» 

430 REM MAM IS LABEL 

440 FOR I5-1S-T7 TO H3 STEP 15 

450 M45SH3+I5-1 

4?0 PLOT 5*I*^. L/1 

500 TO«lOO*H3+M4+4O#IHTiM4/S0> 

510 LABEL >:». 1,5. 1.7.0) 

500 OPLOT “0. 3.-0. 5 

530 LABEL <540. 1 . 2. 1 , 7. 270. 8/1 1 )T6 

540 FORMAT F5.0 

550 tlEMT 15 

5i50 PEM AMIS DONE 

850 GOSOB 330 

860 PRINT "NEN PLOT? PREPARE PLOTTER" 
3?0 PRINT "1«YES" 

800 INPUT R 

900 IF Ral THEN 190 

910 STOP 

900 PEM POWER VS HEIGHT 

93U PEM IHIT. T ARRAY 

940 FOR K«1 TO 128 

950 TCKIeO 

960 NEXT K 

980 FOR 11=1 TO H3 

990 OFFSET 5»I1.0 

1000 REM PLOT ONE MINUTE 

1010 REM FIRST POINT 

1020 C2»-l 

1030 I2.«0 

1040 G0SU8 1900 

1050 Zl-22 

1060 X»<21-B)»H 

1070 Y=5 

1080 GOSUB 1600 

1090 REM REMAIHIHC POINTS 

1100 FOR I2«.l TO 19 

1110 GOSUB 1900 

1120 S=<22-21)/5 

1130 REN S=M INC. BETWEEN PLOT POINTS 

1140 FOR 13=1 TO 5 

1150 V= 1 2*5+ 13 

1160 ;-:r,<Zl-B+I3*S)*H 

1170 GOSUB 1600 

1180 NEXT 13 

1190 21«22 

1200 NEXT 12 

1210 PEN 

12S0 NEXT II 

1230 OFFSET 0.0 

1240 XAXIS 108.H.5O0.5OO+H 

1250 PLOT 50O+<H. 2). loa. 1 

1260 LABEL <#.1.2.1.7.0.8/11) 

1270 CPLOT “2.3- I 

1280 LABEL <#> 1 . ,v.> k . 7. 0. 8/1 1 ) " 10 DB" 
1290 RETURN 
1572 STOP 

1580 REM C2=-l FOR lUST HIDDEN 
1590 REM C2=l FOR Lh'T NOT HIDDEN 
1600 IF <TCY)-5)>I' THEN 1660 
1610 REN NOT HIDDEN 
1620 TrYI=X 
1630 Cl =2 
1640 GOTO 1680 
1650 PEN HIDDEN 
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teso TrV>TrV3-5 
J670 CJ-1 
i680 C3«C1»C2 
1690 C2«t-l>tCl 
1700 PLOT X»Y*C3 
1710 RETURH 

1715 REH HEIGHT FACTOR H DETERtllHES THE 110, OF HORIZONTAL PLOT UHITO 

1716 REM FOR 10 DB. NOTE THAT H0RI2. TICKS ARE 25 PLOT UNITS APART. 
1720 PRINT "INPUT HEIGHT FACTOR" 

1730 INPUT H 

1740 REM BASE VALUE B IS SUBTRACTED FROM ALL DATA BEFORE PLOTTING. 

1741 REM B SHOULD BE NEAR THE MIN, 

1750 PRINT "INPUT BASE VALUE" 

1760 INPUT B 

1770 REM LIMIT VALUE IS THE LIMIT IN DB ABOVE THE MIN, FOR 

1771 REM WHICH THE PLOT WILL BE CLIPPED IF A 1 MIN, SPIKE OCCURS. 
17S0 PRINT "INPUT LIMIT VALUE IN DB" 

1790 INPUT L3 
1800 L3»INTa3) 

1810 RETURN 

1G20 STOP 

1830 H»50 

1840 B«S2+0.08 

1850 L3"A64-0,3 

I860 RETURN 

1870 STOP 

1890 REM LIMITING 

1900 22-DC I2+l» II VlUO 

1910 IF I1-N3 THEN 1940 

1920 IF 22 <« T5 OR DCI2+1 , 1 l + nz-lOO >- T5 THEN 1940 
1930 22»T5 
1940 RETURN 
1550 END 
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